General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



TP 






' 


Final Report 


. 


■ 

, - 














FINAL REPORT 
EE-MSFC- 1900 


CALCULATIONS FOR INTERPRETATION OF 
SOLAR VECTOR MAGNETOGRAPH DATA 


By 

A. R. Dunn, Ph. D. 


May 1975 


Prepared For 

ELECTROMAGNETIC AND SOLID STATE DIVISION 
SPACE SCIENCES LABORATORY 
GEORGE C. MARSHALL SPACE FLIGHT CENTER 
HUNTSVILLE, ALABAMA 


Contract No. NAS8-26376 


Prepared By 

RESEARCH DEPARTMENT 
ELECTRONICS AND ENGINEERING 
TELEDYNE BROWN ENGINEERING 
HUNTSVILLE, ALABAMA 



ABSTRACT 
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ANALYSIS OF SUNSPOT 
SPECTRAL DATA 
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INTRODUCTION 

Observations of the profile of the neutral iron line at 0. 5250216 
micrometer (5, 250.216 A) are described in this section. The reduction 
of the observations to obtain residual intensities, line widths, and the 
locations of the Zeeman components is discussed. The components of 
magnetic field strength along the line of sight are determined and plotted 
in the form of a sunspot map. Finally, the steps to be followed in a 
more complete analysis of the data are outlined. 




SUNSPOT SPECTRUM OBSERVATIONS 

This section describes the analysis of sunspot spectral observa- 
tions made during July 1969 by Dr. M. J. Hagyard with the spectrograph 
of the McMath solar telescope at the Kitt Peak National Observatory. 

(The spectrograph and telescope are described in Reference 1). Although 
both photographic and photoelectric observations were made, the photo- 
graphic data are of primary interest here. 

Because the spectra were intended for Zeeman effect studies, the 
polarization form of the observed intensity had to be known. A partial 
linear polarization is introduced within the telescope in the reflection at 
the primary mirror; this was compensated in the observations by placing 
a glass flat in the light path. The only remaining instrumental effect is 
then a possible phase shift introduced at the primary. Such a phase shift 
would not be detectable in the results before a more advanced stage in 
the analysis is reached. In order to check for residual instrumental 
polarization, photoelectric and photographic observations were made of 
the photosphere, which should be intrinsically unpolarized. Photospheric 
intensities were measured in right and left circularly polarized light and 
in light linearly polarized at 0, 45, and 90 degrees to the entrance slit. 

No dependence of intensity upon polarization form was observed for any 
of the photospheric observations; residual instrumental polarization is, 
therefore, negligible. 

Control over the polarization form of the observed light was 
achieved by placing polarizers in front of the spectrograph slit. The 
slit was always parallel to the east -west direction on the Sun. Immedi- 
ately preceding the slit was a calcite analyzer oriented to pass only those 
light components which were linearly polarized parallel to the slit. 

The calcite analyzer was preceded by a soleil compensator, which can 
be adjusted to form a quarter-wave or a half-wave plate. The form of 
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polarization transmitted is defined by specifying the angle a between the 
compensator's axis and the transmission direction of the calcite analyzer. 
The polarization forms transmitted by the compensator-analyzer com- 
bination are summarized in Table 1; the sequence numbers 1 through 5 
were used to identify the combinations throughout the observations. 

A predisperser was used to eliminate overlapping orders. The 
slit width was 0. 1 millimeter for all observations; the effective slit 
length was determined by the soleil compensator setting. The exposure 
times were either 6 or 12 seconds, and stepwedge calibrations were 
made for both these exposure times. All the spectra discussed in this 
report were photographed on 14 July 1969 between 2:40 and 3:15 p, m. 

M. S. T. The sunspot observed was then located near the center of the 
dis>k (cos 0 = 0. 939). The spot penumbra was approximately circular, 
but the umbra was divided by a light bridge. The position of the umbral 
center was marked on the spectra by a wire placed perpendicular to the 
slit. A second wire was placed parallel to the slit near the red end of 
the observed spectral range for a permanent position marker. 


The spectrograph slit was first placed just south of the sunspot, 
with the length of the slit directed east-west on the Sun. A sequence of 
five exposures was made in the order given in Table 1. The slit was 
then moved 1 millimeter toward the north pole of the solar image, and 
the exposure sequence was repeated. The procedure was repeated, the 
slit being moved north by 1 millimeter each time, until the entire spot 
had been traversed. The solar image formed by the McMath telescope 
is approximately 800 millimeters in diameter, so a displacement of 


1 millimeter in the image corresponds to approximately 1, 740 kilometers 
on the Sun. Twelve positions on the Sun were observed with an expos- 


ure time of 6 seconds. The three positions containing the sunspot umbra 


were also photographed with a 12 -second exposure time. 
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TABLE 1 . OBSERVED POLARIZATION FORMS 


SEQUENCE 

NUMBER 

SOLEIL 

COMPENSATOR 

SETTING 

a (de<j) 

FORM OF 

TRANSMITTED LIGHT 

1 

X/4 plate 

+45 

Left circularly polarized 

2 

X/4 plate 

-45 

Right circularly polarized 

3 

x/2 plate 

0 

Linearly polarized parallel to slit 

4 

X/2 plate 

+22.5 

Linearly polarized at 45 deg to slit* 

5 

x/2 plate 

+45 

Linearly polarized at 90 deg to slit* 


*Recall that a half-wave plate alters the plane of vibration by 2a. 











The photographic spectra were converted to digital form by 
measuring the photographic density as a function of position on the film 
with a digitizing microdensitometer. The films cover approximately 
5 X 10 -4 micrometers (5 A) of the solar spectrum centered at 0. 525 
micrometer (5,250 A); a shorter region of approximately 1. 5X 10" 4 
to 2 X 10“ 4 micrometers (1.5 to 2. 0 A) containing the neutral iron 
line at 0. 5250216 micrometer (5, 250. 216 A) was measured. The 
position scale was converted to dispersion in A/mm by referring to 
the locations of the spectral lines and the wire position marker. The 
microdensitometer sample spacing could then be expressed in terms 
of angstroms on the film; the sampling points were found to be spaced 
a't 5.88 X 10" 7 micrometer (5. 88 mA) intervals. 

The films were sampled in a uniform pattern, which simplifies 
the correlation of the scans with their proper positions in the sunspot. 
The pattern is shown in Figure 1. 


The spectral dispersion was along the north- south direction, 
so each microdensitometer scan was made parallel to the dispersion 

i 

and had a height equal to the microdensitometer slit height of 1 milli- 
meter. The width of spectrum sampled at each scan step was deter- 
mined by the microdensitometer slit width of 30 micrometers. At the 
dispersion of the film, this represents a sample width of 4. 20 X 10" 
micrometers (4.20 mA). 


v-7 


The center of the umbra is marked on the film by the image of 
the wire placed perpendicular to the slit. Microdensitometer scans 
were made beginning at the wire and working eastward, the final scan 
being entirely in the photosphere. These scans were labelled a through 
f. Scans were then made westward from the wire into the penumbra; 
these were labelled g through i. For those exposures which did not 
contain the entire sunspot, the same pattern was followed, but the 
reference line was the centerline of the film instead of the center of 
the umbra. 
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FIGURE 1. MICRODENSITOMETRY PATTERN FOR SUNSPOT OBSERVATIONS 







The digital intensity data are labelled with three symbols which 
describe the position of the observed region and the polarization form 
of the recorded intensity. For example, in the label X-2-c, the Roman 
numeral X designates the north- south location of the region observed, 
in terms of millimeters on the primary solar image; 2 is the polariza- 
tion form, in this ease right-hand circular polarization; c is the east- 
west location of the measured strip, in terms of millimeters on the 
film. There are 564 strips in the 15 exposures made on 14 July 1969, 
but not all of these contain sunspot spectra. 

The microdensitometer output represented measured film density 
in the sample areas of the film. The digital output data were recorded 

r 

on magnetic tape and were converted to intensities by use of the charac- 
teristic curves derived from the stepwedge calibration exposures. The 
magnetic tapes thus contain intensities uniformly spaced 5. 88 X 10“ 7 
micrometers (5. 88 mA) apart, defining a region approximately 1 . 5 X 
10" 4 to 2 X 10" 4 micrometers (1. 5 to 2. A) long. The intensity values 
are in arbitrary units, in the sense that the continuum level shown on 
the tape records is usually considerably less than 1. 00. The intensity 
values are in correct proportion to one another but should be rescaled 
to bring the continuum up to 1. 00. These digital data were the values 
used in the analysis to be described in the subsection entitled "Line 
Profile Data Reduction Program". 
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REDUCTION OF THORIUM CALIBRATION WAVELENGTHS 
TO STANDARD CONDITIONS 

To obtain a laboratory determination of the wavelengths of the 
solar spectrum lines, a thorium emission spectrum and a solar disk 
spectrum were exposed in parallel rows on the same film. The 
observed solar wavelengths can then be determined by measuring their 
distances from thorium standard reference lines (reference standards 
are discussed in References 2 and 3; precise determinations of addi- 
tional thorium wavelengths are given in References 4, 5, and 6). This 
procedure refers the solar lines to the wavelengths for the thorium 
lines at the time of observation; the observed wavelength actually 
depends slightly upon the conditions pf observation. 


Because the index of refraction of air depends upon its tempera- 
ture and pressure, and upon the wavelength of radiation being studied, 
measured wavelengths are reduced to the values they would have in 
some chosen standard conditions. The usual standard conditions are 
dry air at temperature 15° C, pressure 760 millimeters of mercury; 
or vacuum. Air values are used more often as standards for wave- 
lengths in the visible range. 


In 1952, a meeting of the Joint Commission for Spectroscopy, 
sponsored by the International Astronomical Union, the International 
Union of Pure and Applied Physics, and the International Council of 
Scientific Unions, recommended adoption of Edlen's formula for con- 
verting wavelengths in standard air to wavelengths in vacuum (Ref. 7). 
In the same paper, Edlen also gives a simple formula for converting 
measured wavelengths to values for standard air. This simple formula 
can be used when the conditions of measurement are not very far from 
the standard conditions. An additional correction for the presence of 
water vapor should be made in high-precision work. 
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A short FORTRAN program (Figure 2) has been written to correct 
wavelengths to standard air using Edlen's formula. The partial pressure 
of water vapor was not measured with the wavelengths to be corrected, 
so a moisture correction could not be included. The program uses the 
following version of Edlen's formula: 



AX. i 


X 2 \ /0. 0013882 p 
\i/\l + 0.00367 t 



where 


o 

X2 


X-2 


Xl ! 

AXi , AX2 

P 

t 


unknown wavelength in standard air 
measured value to be corrected 
reference wavelength for standard air 
vacuum corrections to \i and X 2 
air pressure, mm Hg 
temperature, °C. 


The vacuum corrections , AX 2 are tabulated in Table 3 of 
Reference 8. 

The program has been tested on the two examples of wavelength 
corrections given by Babcock (Ref. 9). The results of the test are 
summarized in Table 2. 


TABLE 2. RESULTS OF TESTING EDLEN CORRECTION PROGRAM 


p (mm) 

t (°C) 

*1 (u ) 

\ _2 (m ) 

BABCOCK (v ) 

DUNN (y) 

620 

20 

0.4500 

0.300 

-6.4 x 10- 7 

-6.6 x 10- 7 

720 

25 

0.400 

0.800 

+5.2 x 10“ 7 

+5.2 x 10“ 7 
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REDUCTION OF WAVELENGTHS TO STANDARD AIR 10/22/70 

REFS, •_ INTRO. TO TABLE OF WAVENUMBERS 2000A TO 7000 A > N.B.S. 
MONOGRAPH 3* VOL. 1 (GET VACUUM CORRECTIONS FROM 
TABLE 3 OF THIS REF.)* EDLEN * JOSA VOL. 43* PG. 339* 1953* 
BABCOCK* AP. J.* VOL. Ill* PG. 60 * 1950. 

EDLEN FORMULA 

MUST TREAT WAVELENGTHS IN 2 PARTS TO MINIMIZE LOSS OF 
SIGNIFICANCE. CHARACTERISTIC WILL BE READ AS INTEGER 
VARIABLE* MANTISSA AS REAL VARIABLE. LW1C* W1M ARE REFERENCE 
WAVELENGTH* D1 IS ITS VACUUM CORRECTION. LW2C » W2M. D2 ARE 
SAME QUANTITIES FOR WAVELENGTH TO BE CORRECTED. 

D2 IS THE TABULATED VALUE FOR THE MEASURED WAVELENGTH. D2* D1 
AND REF. WAVELENGTH ARE VALUES FOR STANDARD AIR. 

NEGLECT MOISTURE CORRECTION. 

P IS PRESSURE IN MM HG * T IS TEMPERATURE IN DEGREES C. 

K P U = 2 
KPR =3 

READ ( KPU ♦ 1) LW1C * W1M* D1 

1 FORMAT ( 16* 2F10.0) 

WRITE (KPR* 3) 

3 FORMAT (1H *3 7H WAVE LENGTHS CORRECTED TO STANDARD AIR) 

WRITE (KPR* 7) 

7 FORMAT (1H * 37HADD CORRECTION TO MEASURED WAVELENGTH) 

WRITE (KPR* y) LW1C* W1M 

5 FORMAT (1H * 2 1HREFERENCE WAVELENGTH® ♦ 16* 1H+* F10.6) 

WRITE (KPR* 4) 

4 FORMAT (1H * 19HME ASURE D WAVELENGTH* 6X » 1 5HCORRECT I ON IN A) 

READ (KPU* 2) P. T 

2 FORMAT ( 2F10.0) 

10 READ (KPU* 1) LW2C * W2M * D2 
LOOK FOR LAST CARD (BLANK) 

IF (LW2C-0) 6* 60* 6 

6 W2C - LW2C 
W1C* LW1C 

A* W2C / ( W1C + W 1 M ) 

B= W2M/( W1C + W 1 M ) 

FAC* ( 1.3882E-3*P/( 1.0 + T*3 .67E-3 )>-l .0 
CORR= (D2-D1*(A + B))*FAC 
WRITE (KPR* 8) LW2C* W2M* CORR 

8 FORMAT (1H ♦ 3X ♦ I 6 * 1H+ *F10.6 * 8X * F10.6) 

GO TQ^.l 0 

60 CALL EXIT 
END 

FIGURE 2. EDLEN CORRECTION TO MEASURED WAVELENGTHS 
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In practice, the corrections derived from this program will not 
be required. Unless the measurement conditions deviate grossly from 
standard, the wavelength corrections will be very small for wavelengths 
which lie within approximately 5 X 10" 3 micrometers (50 A) of the line 
used for a wavelength reference. The observed thorium calibration 
spectrum is very short, covering the range 0.5248 to 0. 5256 micro- 
meter (5, 248 to 5, 256 A* approximately. There is no problem here 
of measuring wavelengths from a reference standard hundreds or thou- 
sands of angstroms away. The variation in wavelength caused by labora- 
tory air will be negligible over the short observed wavelength range. 




LINE PROFILE DATA REDUCTION PROGRAM 


The sunspot spectra are to be used to determine sunspot magnetic 
field strengths from the profiles of the magnetically sensitive neutral 
iron line at 0. 5250216 micrometer (5, 250. 216 A). The digital intensity 
data described in the first subsection require additional reduction to 
prepare the profiles for determination of magnetic field strengths. A 
FORTRAN program (Figure 3) has been written to perform the reduction. 

The intensity data are used here without correction for the 
effect of stray light. Stray light strongly distorts measured umbral 
continuum intensities, but the distorting effect is much smaller in the 
case of residual intensity in a spectral line which is present in both the 
sunspot and the photosphere. The 0.5250216 micrometer (5,250. 216 A) 
line is, in fact, stronger in spot umbrae than in the photosphere. The 
observed sunspot profiles without stray light correction are generally 
sitronger than the photospheric 0. 5250216-micrometer (5,250. 216 A) 
line in the Utrecht Atlas; and it was thought that in this case the distor- 
tion of the profiles by penumbral and photospheric stray light was 
smaller than the distortion by such other causes as blending of Zeeman 
components, absorption by telluric water vapor, and blending with the 
neighboring faint molecular line. 

The digital intensity values appear on the magnetic tape in order 
of decreasing wavelength; the values are spaced 5. 88 X 10“ 7 micrometers 
(5. 88 mA) apart. The spectral range runs from the marker wire near 
0. 5254 micrometer (5,254 A), to approximately 0. 5249 micrometer 
(5,249 A). The portion of spectrum between 0. 5252 and 0. 5253 micro- 
meter (5,252 and 5,253 A) reaches the true continuum in the photosphere, 
according to the Utrecht Atlas. Examination of the sunspot spectral data 
shows that the measured intensities are generally largest in the 9. 5252 
to 0. 5253- micrometer (5, 252 to 5, 253 A) interval, so these intensities 
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PROFILE PARAMETERS 12/15/70 ARD 

7094 VERSION 1/5/71 

FIND MINIMUM VALUES* HALFWIDTH* AND EQUIVALENT WIDTH OF 5250.2 
PROFILE. 

Y * S ARE UNSCALED INTENSITIES* CONT IS CONTINUUM LEVEL FOR SCAN. 

RLAM IS RESIDUAL INTENSITY* W IS EQUIVALENT WIDTH. LAST Y 
VALUE WILL BE FOLLOWED BY Y= 2.0 TC SIGNAL END OF DATA. ID IS 
IDENTIFICATION NUMBER FOR SCAN* IBC IS 'BROADENING CODE' 

DESCRIBING EXPECTED FORM OF ZEEMAN BROADENING. 121 IS LEFT OF 
PROFILE CENTER. IT IS THE RIGHT HAND POINT OF ITS INTERVAL. 

12? IS THE LEFT HAND POINT OF THE SAME INTERVAL. Ill AND 112 
ARE RESPECTIVELY THE LEFT AND RIGHT HAND POINTS OF AN INTERVAL TO 
THE RIGHT OF THE PROFILE CENTER. 

DIMENSION Y ( 1 0 ) * RLAM(150>* RMINC6) 

DIMENSION XMINC6) 

DIMENSION NUN ( 5 ) * FT ( 3 ) * FNNT ( 3 ) • FOT < 2 ) 

DIMENSION FAT ( 2 ) * MALL ( 2 ) 

DIMENSION ILOSE ( 9 ) 

DATA FOT / 8H ( I 1*2A6> / 
data FNNT / 18H ( 2 A6 • I 3 * I 5 » IX *9A6 ) / 

DATA FAT/5H15A6)/ 

DATA FT/15H( I 1 .F6.4.9F7.4) / 

KPU = 5 
KPR= 6 

WRITE (KPR* 15) 

WRITE ( KPR * 16) 

15 FORMAT ( 49H POSITIONS ARE IN TERMS OF DATA SPACING INTERVALS) 

16 FORMAT (26H AN INTERVAL IS ABOUT 5 MA ) 

C AN INTERVAL IS 5.88 MA FOR DATA OF 14 JULY 1969 

300 CALL REDTPD ( 8 *FOT * I ER • 1 * ICK* 2* MAL'L) 

IF HER .EQ. 2 ) GO TO 6 7 
IF (ICK . NE. 1) GO TO 300 
C LOOK FOR DATE. (FIRST RECORD IN SCAN) 

68 READ (KPU* 1) ID. IBC,. IRR* CONT 
1 FORMAT (315. F10.0) 

C IBC IS NUMBER OF RESOLVED AND UNRESOLVED MINIMA OBTAINED FROM VlSLAu 

C INSPECTION OF PROFILE PLOTS. IT IS USED ONLY AS A MEANS OF i. 

C CHECKING THE PROGRAM. WHICH DETERMINES THE NUMBER* TYPES. j 

C AlkD POSTIONS OF MINIMA IN A DIFFERENT WAY. IBC= (NO. OF RESOLVED 

C MIN. ) * 1 0 + NO. OF UNRESOLVED MIN. 

IF (ID) 67* 67* 301 

C NOW READ 2ND RECORD AND CHECK TO MAKE SURE YOU HAVE THE RIGhT 

C SCAN 


FIGURE 3. LISTING OF LINE PROFILE REDUCTION PROGRAM j 


17 


nnn r\ n n 


•*01 CALL REDTPD (8* FNNT * I ER » 2* NUN * if I ROLL t if I PAP* 9* I LOSE) 
If UROLL *EQ« IRR) GO TO 30 
18 CALL REDTPD 18* FOT • I ER * 1* ICK* ?* MALL) 

IF ( I ER • EQ* 2 ) GO TO 67 
IF ( ICK .NE. 1 ) GO TO 18 

CALL REDTPD i 8 *FNNT * I ER *2 *NlJN * 1 * I ROLL * 1 * I PAP *9 * ILOSE) 

IF (I ROLL *NE. IRR) GO TO 18 

IF IRR NOT EQUAL TO I ROLL# SAVE IRR AND CONT AND SKIP TO NEXT 
SCAN 

30 CALL REDTPD (8# FAT# I ER t St NUN) 

66 J = 1 

J IS USED TO COUNT STORE VALUES 
WRITE ( KPR t 12) ID* IROLL 
1? FORMAT (1H * 8HID* NO.** 15* 7H IROLL** 15) 

YT * CONT * *475 

5 CALL REDTPD (8* FT* I ER * 1* ICK* 10* Y I 1 ) ) 

? FORMAT (10F7.4) 

DO 7 I* 1* 10 
IF (Y(I )-YT) 4* 7* 7 
4 IF ( Y I I ) ) 300# 300* 9 
9 STORE = Y ( I ) 

WRITE (KPR* 74) STORE 
74 FORMAT < 1H • F10.5) 

73 IF ( I -10) 71# 70* 70 

70 CALL REDTPD (8* FT* I ER • 1* ICK* 10* YU)) 

I* 1 

GO TO 72 

71 I* I + 1 

72 IF (Y(l)-YT) 73* 73* 104 
104 IF (J-3) 17* 8* 6 

17 J ,* J + 1 

7 CONTINUE 
GO TO 5 

8 K s 1 

DO 119 KK = 1# 15 

WRITE THE 150 POINTS AFTER THE THIRD ‘STORE* INTO RLAM ARRAY 
BUT START WITH 1ST DATA POINT OF NEXT GROUP OF 10* NOT NEXT 
DATA POINT 

CALL REDTPD (8* FT* I ER * 1* ICK* 10# Y ( 1-1 > 

DO 6 I* 1* 10 
RLAM ( K ) = Y ( I ) /CONT 
K= K + 1 

6 CONTINUE 
119 CONTINUE 

11 CALL REDTPD (8* FOT* IER* 1* ICK* 2* MALL) 

IF (IER .EQ. 2) GO TO 67 
IF (ICK *NE. 1) GO TO 11 


FIGURE 3 - Continued 
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C FIND MINIMA WITHIN 5250.2 

13 ICHK* 0 

C I CHK IS VARIABLE USED TO COUNT MINIMA 

no • io 

n= l 

DO 250 KK = It 15 

WRITE (KPR* 251) (RLAM(I)t 1“ lit 110) 

251 FORMAT ( 1H * 10F7.4) 

11= II + 10 

iio = no + io 

250 CONTINUE 
I * 40 
NR= 0 
NUR = 0 

81-DELM# RLAMU )-RLAM( 1-1 > 

DEL= RL AM ( 1 + 1 ) - RL AM ( I ) 
r DELP= RLAMU+2) - RLAMII+1) 

IF (DELM) 83* 82* 83 

82 1= I + 1 
I GO TO 81 

83 IF ( DELP ) 84* 88 * 84 
88 IF (DEL) 89* 90* 89 

90 IF (RLAM( I+3)-RLAM( 1+2) ) 89* 89* 91 

91 I CHK* ICHK + 1 
RMIN ( I CHK ) = RLAMt I -*-l > 

XMIN(ICHK) =1+1 
1 = 1 + 1 
NUR= NUR + 1 

WRITE (KPR* 155) ICHK* XMINUCHK)* RM I N ( I CHK ) * NR. NUR 
IN THIS CASE* WANT TO ARRIVE FINALLY AT STATEMENT NO. 82 
WITH I« ORIGINAL I + 2 • SO INCREASE I TWICE 
155 FORMAT (10H MIN. NO. • I 5 * 5X * 2H I * *F6 .2 *5X * 5HRLAM- * F8 .4 » 2 I 7 ) 
GO TO 89 

84 IF (DEL) 85* 92* 85 

85 IF ( ABS ( DELM + DLL ) - (ABS (DELM) + ABS (DEL))) 86* 87* 87 

92 IF (DELM) 93* 89* 89 

93 IF (DELP) 89* 89* 94 

94 ICHK* ICHK + 1 
RM I N ( ICHK) = RLAM ( I ) 

XS= I 

XMINUCHK) = XS + 0.500 
NR= NR + 1 

WRITE (KPR* 155) ICHK* XMIN(ICHK)* RMIN ( ICHK ) * NR. NUR 
GO TO 89 

86 IF (DELM) 96* 89* 89 

96 IF (DEL) 89* 89. 97 

97 ICHK- ICHK + 1 
I SV * I 


FIGURE 3 - Continued 

ORIGINAL PAGE IS 

OF POOR QUALITY lq 


IS* 


non 


NR = NR + 1 

RMINI ICHK )= RLAM ( I ) 

XMIN(ICHK) = I 

WRITE ( KPR* 155) ICHK* XM.I N ( I CHK ) t RMINIICHK)* NR* NUR 
GO TO 89 

87 IF (ABS(DEL) - ABS (DELP)) 98* 89. 89 

98 IF ( ABS t DEL ) - ABSCDELMM 99* 89* 89 

99 ICHK * ICHK + 1 
RMIN(ICHK) = RLAM ( I ) 

XMIN(ICHK) * I 

NUR = NUR + 1 

WRITE ( KPR • 155) ICHK* XMINUCHK)* RM I N ( I CHK ) » NR* NUR 
89 IF (1-65) 82* 100* 100 

100 IF (RLAM(I) - *800) 82. 82. 101 

101 I TRY* 1 0*NR *► NUR 

IF (IBC-ITRY) 103* 102* 103 
103 WRITE (KPR. 105) NR* NUR. ITRY. I BC 
105 FORMAT ( 4H NR* * I 5 • 3X *4HNUR= * 1 5 * 3Y * 5H I TRY* * I 5 ♦ 3X * 
116HBROADENING CODE* * 15) 

WRITE (KPR* 107) 

107 FORMAT (25H CHECK DATA FOR THIS SCAN) 

102 IF (ICHK-3) 19* 20. 20 

19 ICHK * ICHK + 1 
RMIN(ICHK) = 1.00 
GO TO 102 

20 IF (RMIN(l) - RM I N ( 2 ) ) 21* 21* 22 

21 ABSMN* RMIN(l) 

GO TO 23 

22 ABSMN* RM I N ( 2 ) 

23 DO 924 I* 3. 6 

IF ( RM I N ( D-ABSMN) 57* 57* 924 
5 7 ABSMN* RM I N ( i ) 

924 CONTINUE 

WRITE (KPR* 128) ABSMN 
128 FORMAT ( 1H * F10.4) 

C START HALFWIDTH CALCULATION 

24 HMIN* 0 • 50* ( 1*0 ♦ ABSMN) 

11= ISV 

ISV IS LOCATION OF FIRST RESOLVED MIN. THERE MUST ALWAYS BE 
AT LEAST ONE OF THEM. 

WORK TO RIGHT OF PROFILE CENTER FIRST . 

25 test* hmin - rlami i i ) 

IF (TEST) 59* 58* 58 

58 11= 11 + 1 
GO TO 25 

59 112* 1 1 
111* II - 1 

C NOW WORK TO LEFT OF PROFILE CENTER 

12* ISV 

FIGURE 3 - Continued 
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26 TEST* HMIN - RLAMU2) 

IF (TEST) 48 • 49 . 49 

49 12® 12 - 1 
GO TO 26 
48 122® 12 

j 121* 12 + 1 

FRAC* (HMIN-RLAM( 121) > / ( RLAM < 1 22 ) -RLAM ( 121) ) 

LINEAR INTERPOLATION TO GET POSITIONS OF HALFWIDTH POINTS 
THIS SHOULD BE OK BECAUSE PROFILE SHOULD BE NEARLY LINEAR IN 
THIS RANGE 
X I 2 1 * 121 
POS2* X I 2 1 - FRAC 

FRAC* (HMIN-RLAMl 111) )/(RLAM( 1 12) -RLAM ( 111) ) 

XIII* 111 
POS1* XI 11 + FRAC 
FWHM* POS1 - POS2 
WRITE ( KPR » 60) P052. POSlt FWHM 
60 FORMAT ( 1H * 1 9HHALFW I DTH POINTS AT. F7.3.5H AND.F7.3t 
17H FWHM*. F7.3) 

C START EQUIVALENT WIDTH CALCULATION 

I® 122 

C LOOK FOR ENDS OF PROFILE 

27 X* RLAMU ) -RLAM (1*1 ) 

IF (X) 46* 47. 47 

46 I * I - 1 
GO TO 27 

47 ISTRT* I 

WRITE (KPR. 147) ISTRT 

147 FORMAT (1H * ISTRT* . 15) 

1*112 

? 8 X * RLAM ( 1+1 ) - RLAM ( I ) 

IF (X) 45. 45. 44 

44 I* I + 1 
GO TO 28 

45 ISTP* I 

WRITE (KPR. 148) ISTP 

148 FORMAT ( 1H ♦ 5HISTP*. 15) 

NI* ( ISTP-ISTRT + 1 >/2 
XI* ISTP - ISTRI ♦ 1 

X I 2 = X I / 2 • 

XI* NI 

IF ( X I —X I 2 ) 2 9. 43* 29 
43 ISTP* ISTP - 1 
29 SUM* 0* 

DO 41 I* It 150. 2 
C WILL NEVER GET TO I * 150 

IS * ISTRT ♦ I 
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1 SUM* RLAM(IS>*4* + RLAM(IS + 1>*2. + SUM 
IF ( IS -< ISTF - 3 )> 41* 62# 41 

41 CONTINUE 

62 $* RLAMUSTRT) + RLAM(ISTP) 4 4.*RLAMiISTt> - 1) 

SUM* (5 + SUM ) / 3 • 

DLAM* ISTP - ISTRT 
W= DLAM - SUM 
WRITE <KPR* 65) W 

65 FORMAT < 3H W* *5X *F10«5 *5X *22H IN DATA SPACING UNITS) 
GO TO 68 

67 WRITE ( KPR • 10) 

10 FORMAT <9H END FILE) 

STOP 

END 


FIGURE 3 - Concluded 

ORIGINAL PAGE IS 

OF POOR QUALITY . 


were adopted as the continuum values for the sunspot scans, and the 

i 

remaining intensities were scaled to them. The continuum intensity 
was determined for each scan by inspection of the magnetic tape print- 
out, and that value was read as the variable "CONT” for the line profile 
reduction program. 

Starting at the red end of the region, and including the image of 
the wire, the measured intensity falls below 50 percent of the continuum 
intensity three times before the 0.52 50216-micrometer (5, ".50.216 A) 
line is reached. This fact was used to isolate the 0. 5250216-micrometer 
(*>, 250.216 A) profile from the rest of the spectrum. The minima deeper 
than some fraction of CONT were counted, and the third minimum was 
known to be X.5250. 6, the strong line immediately adjacent to \5250.2. 

(The fractional value used was usually 0.475. However, because the 
relative ^tensities of the spectral lines are different on different scans, 
this value was sometimes adjusted to avoid counting unwanted lines or 
overlooking one of the desired minima. It was always possible to find 
a fraction which would select the three desired minima. ) 

When the third minimum has been found, the program tests the 
succeeding intensity values until it finds the next local intensity maximum, 
which represents the dividing line between the X5250. 6 profile and the 
X.5250.2 profile. (The intensity usually does not quite return to the 
continuum level between the two profiles. ) The next spectral line, 
then, will be X5250.216. An array of 150 values is set aside for that 
profile. One hundred fifty intensities following the local maximum are 
read into this array, called RLAM, and the rest of the unsealed digital 
values are not used. The intensities are scaled by dividing them by 
CONT before they are stored; they are stored as residual intensities. 
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The program next locates the intensity minima within the 
0. 52 50216-micrometer (5, 250.216 A) profile. * These minimum value 
locations will be used later as the positions of the Zeeman components. 
The program begins searching at the fortieth point in the profile. In a 

i 

few cases, unresolved minima appeared before the fortieth point, so 
the starting index should probably be changed. The presence of such 
neglected minima can always be detected because the program prints 
the 150 RLAM values for each profile. 

The program searches for two kinds of minima, called "resolved" 
and "unresolved". A "resolved" local minimum is defined as an intensity 
value which lies between two greater intensities. The intensity value 
immediately adjacent on one side may be equal to the minimum value. 

For example, in the sequence 0. 52 10, 0. 5009, 0.4841, 0.4995, 0. 5196, 
the value 0.4841 would be selected as a resolved minimum. Its position 
would be given by its index: in this case, 3. In the sequence 0. 5210, 
0.5009, 0.4841, 0.4841, 0.4995, the resolved minimum would again be 
the lowest intensity value, 0.4841, but the position would be taken as 
halfway between the two intensity points with that value, or 3. 5. An 
"unresolved" local minimum is defined as the left-hand point of an 
interval which is bounded by two intervals having slopes of the same 
sign but greater magnitude. Because the digital data are smooth, the 
use of minima found over a region only four or five points long is 
successful; noise fluctuations are rarely detectable on this scale. 

As a check, the number of resolved and unresolved minima is 
also estimated visually from plots of the profile scans, and a code 


*A11 positions found by this program are the array indices of the corres- 
ponding intensity values. They may be converted to wavelength values 
by using the fact that adjacent intensity values are separated by 5. 88 
X lO’ 7 micrometers (5.88 mA). 
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number (IBC) equal to ten times the number of resolved minima plus 
the number of unresolved minima is read for each scan. This code 
number is checked against the number of minima found by the program, 
and a warning message is printed if they do not agree. Lack of agree- 
ment is not necessarily a sign or error, since the program's selection 
criteria are not exactly the same as those of the eye. However, the 
warning message provides a screening mechanism which selects the 
scans that are most likely to contain errors. The data can then be 
checked to ensure that the correct scan is being read and that the 
0.,52 50216-micrometer (5, 250.216 A) line has been correctly isolated. 

The program next determines the full width of the 0. 5250216- 
micrometer (5, 250.216 A) profile at half the minimum intensity, where 
the minimum intensity, I m in» is taken to be the intensity of the deepest 
minimum in the profile. The intensities I^alf at half minimum are then 
equal to (1.00 + I m i n )/2. The program searches outward on both sides 
of the profile from the deepest minimum until it finds the intervals 
containing the intensities Ihalf an< ^ interpolates, if necessary, to find 
the positions of the Ih a if values. 

Finally, the program defines the ends of the profile as the 
positions near the continuum at which intensity reaches local maximum 
values. The indices corresponding to these local maxima are printed. 
The equivalent width W\ is then obtained by direct integration of the 
profile between the end points, using the defining equation 



where 


\ - wavelength 

\i and X 2 - ends of the profile 

I (X) - residual intensity as a function of wavelength. 

The card input and program output are summarized in Tables 3 

and 4. 


TABLE 3. LINE PROFILE REDUCTION PROGRAM INPUT 


VARIABLE NAME 

DESCRIPTION 

ID 

Identification number for program output. Each 
scan has a unique number; they range from 1 to 
564. 

IBC 

Code number for visual estimate of number of 
minima within profile. 

IRR 

; 

Identification number of the scan on the magnetic 
tape. Two scans on different tapes may have the 
same IRR. 

CONT 

Continuum level obtained from the magnetic tape 
printout. 


One card for each scan: ID, IBC, IRR, CONT Format (315, F10.0) 
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TABLE 4. LINE PROFILE REDUCTION PROGRAM OUTPUT 


T) ID, IRR for identification of scan 

2) Intensities of the three minima preceding x5, 250. 21 6 

3) 150-value array containing the profile of A5, 250. 216 

4) For each minimum found, a line is printed containing its position 
and intensity, and a two-digit code defined like IBC, giving a 
running total of the number of resolved and unresolved minima 
found in the profile. 

5) Total number of resolved and unresolved minima, the two-digit 
minimum-count code (ITRY) for the entire profile, and the IBC 
value read, for comparison. 

6) Warning message if ITRY is not equal to IBC. 

7) Intensity of the deepest minimum in the profile. 

8) Locations of the I half points and the full width at half minimum 
in the data spacing units. 

9) and 10) Ends of the profile. 

11) Equivalent width, in data spacing units. 
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MAGNETIC FIELD STRENGTH DETERMINATIONS 
FROM THE 5,250.216 A LINE 

METHODS OF DETERMINING FIELD STRENGTHS 

Determination of sunspot magnetic field strengths from the 
Zeeman splitting of spectral lines is discussed in detail in several 
textbooks (see, for example, Refs. 10 and 11). Reference 10 also 
contains a discussion of Unno's theory of spectral line formation in a 
magnetic field. Unno's technique was one of the first to permit calcula- 
tion of absorption Zeeman profiles in a stellar atmosphere. It has been 
used by a number of investigators and is still used occasionally. It 
assumes a Milne- Eddington atmosphere; profiles which are formed in 
pure absorption and LTE; a uniform magnetic field; and a ratio of line 
to continuous absorption which is constant with depth. 

Two more recent approaches are those of Beckers (Refs. 12 and 

i 

14) and of Moe (Ref. 13). Reference 13 also contains a summary of 
other methods. [A discussion by Staude (Ref. 15) of notational differences 
in recent papers is helpful in interpreting the literature. ] Moe's method 
is probably the most popular one now. It is simple to program and per- 
mits calculation of profiles for any desired atmospheric model. (Refs. 

16 and 17 contain results obtained by Moe's method.) 

Moe's method consists of a different procedure for solving Unno's 
transfer equations for the Stokes parameters. Again, the line profile is 
assumed to be formed in LTE, and the magnetic field is assumed homo- 
geneous over the region of line formation. The requirement of depth- 
independence of the absorption coefficients, however, is relaxed. The 
ratio of line to continuous absorption is represented as a product of a 
depth- dependent, wavelength-independent factor and a depth-independent, 
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wavelength- dependent factor. The emergent intensities of the Stokes 
parameters I, Q, and V can then be written analytically and can be 
evaluated by numerical integration for a given model atmosphere. 
Numerical calculations are practical only for lines formed in pure 
absorption, although Moe indicates the form of the solution for scattering 
lines. 

THE IRON LINE AT 5, 250.216 A 

The iron line at 0. 5250216 micrometer (5, 250.216 A) is often 
us^d for magnetograph field determinations. It has a strong triplet 
splitting and, in general, appears to conform with von Kluber 's criteria 
for suitable lines for magnetic field measurements (Ref. 18). 

Upon closer acquaintance, however, the line has revealed some 
unfavorable characteristics. It is a strong line with a low excitation 
potential and a large scattering component (Ref. 19). Detailed analysis 
is further complicated by the fact that the upper level of the line is 
split. Because the line is strong, it tends to saturate easily; this 
property has been suggested as the reason the line sometimes appears 
to have an "anomalous" Zeeman it - component (Refs. 20 and 21). The 
line is extremely temperature- sensitive, becoming stronger in regions 
of lower temperature (Ref. 22). Moreover, Moe has recently discovered 
a weak molecular line in the red wing of the profile, so the 0. 5250216- 
micrometer (5, 250.216 A) line usually observed is actually a blend. 

Some of the properties of the line are given in Table 5. 


N 

V 
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TABLE 5. PROPERTIES OF FE I A5, 250. 216 


Transition 

t 


i Lower and Upper 
Excitation Potentials 


Zeeman Pattern 


a 5 D 0 - z 7 Dj Ref. 10, p. 179 

0.121 eV, 2.471 eV Ref. 10, p. 179 



Ref. 10, p. 179 


j Lande g Value 3 Ref. 23, p. 382 

i • . : 


fLog g/ -4.46 Refs. 24 and 25 

! (g is statistical weight 
| of lower level ; / is 
oscillator strength for 
.1 transition) 


*Photospheric Fe 
Abundance 



*The proper choice of abundance and /-values of iron is not yet 
established; other published values are as justifiable as those 
given here. (See, for example, Refs. 27, 28, and 29.) 
Fortunately, the g/-abundance product, which is the quantity 
needed for line profile calculations, is better known than 
either of its constituent factors. 
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PRELIMINARY DETERMINATION OF SUNSPOT MAGNETIC FIELD 


For simple triplet splitting, the magnetic field strength along 
the line of sight can be derived immediately from 


AX h 


eH 

4rrm e c 2 



- separation of the it and <r Zeeman components, cm 

- undisturbed wavelength, cm 

- field strength, gauss 

- electron charge, e. s.u. 

- electron mass, g 

- speed of light, cm sec" 1 . 

A simple FORTRAN program has been written to calculate field 
strengths from the triplet splittings obtained from the spectral data. 

The results are field strength components, in gauss, along the line of 

: ■ 

sight. This solution gives no information about the field orientation 

i 

and is limited to determination of field strengths large enough to pro- 
duce resolvable Zeeman components. The field strengths will be 
required for calculating line profiles in the spot; they can also be used 
to construct a crude sunspot map of the projected field vectors. 

The program is given in Figure 4; the results, in the form of a 
sixnspot map, are given in Figure 5. Figure 5 is not drawn to scale. 

The length of a cell in the East-West direction represents approximately 
1, 750 kilometers on the Sun. The length in the North- South direction 


where 


A\ h 

\> 

H 


m e 

c 
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o n n n n r> 


C CALCULATE MAGNETIC FIELD STRENGTHS FOR TRIPLET SPLITTING 

C READ POSITIONS OF MINIMA ND ASSUME THESE ARE CENTERS OF COMPONENTS 

KPU* 2 
KPR = 3 

RFAD ( KPU * 1) W* G 

1 FORMAT (F10.0* ^5.0) 

C WAVELENGTH OF UNSHIFTED LINE IN ANGSTROMS* LANDE G VALUE 

C CONVERT W TO CM 

W = W*1 • OF-8 j 

2 READ (KPU* 3) ID* XL* XC* XR 

3 FORMAT ( 15* 5X. 3F10.0 ) 

READ POSITIONS OF COMPONENTS IN DATA SPACING UNITS 
IF A COMPONENT DOtS NOT APPEAR* LEAVE IT OUT 
IF ALL 3 COMPONENTS ARE PRESENT* ASSUME ZEEMAN SPLITTING IS 
AVERAGE OF THE 2 VALUES THUS DETERMINED. 

THIS PROGRAM ASSUMES AT LEAST 2 MINIMUM VALUES ARE READ' FOR 
EVERY LINE 
NX L = XL 

IF (NXL-0 ) 2# 5* 6 

5 DLH* ( XR-XC)*5.8BE-11 

C OLAMBOAH IN CM 

GO TO 10 I 

6 NXC= XC ■ : 

IF (NXC-O) 2* 7* 8 

7 DLH= ( XR-XL)*2.94E-11 j 

GO TO 10 | 

8 NXRsXR 

IF (NXR-0) 2* 9* 11 

? DLH= < XC-XL)*5.88E-11 j 

GO TO 10 

11 DLH* ( ( XC-XD-MXR-XC) )*2.94E-11 

10 H= DLH/ (4.67t-5*G*W*W) j 

WRITE ( KPR • 12) ID* H I 

12 FORMAT (4H I D* * I 3 * 1QX • 2HH* * F 1 2 • 5 * 7H GAUSS) 

GO TO 2 | 

13 CALL EXIT | 

END 


FIGURE 4. FIELD STRENGTH PROGRAM 
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is approximately 175 kilometers on the Sun, and the separation of 
observations in the North-South direction is approximately 1500 kilom- 
eters. The Roman numerals along the left side of the figure arid the 
lower case letters along the top are the labels of the observations cor- 
responding to the plotted locations. 

DETAILED MAGNETIC FIELD DETERMINATION 

The complete determination of magnetic field strength and 
inclination can be summarized in a diagram (Figure 6). The three 

t 

operations in the left-hand column (observations, determination of 
Zeeman splitting, and calculations of field strength components along 
the line of sight) have been completed. The remainder of the diagram 
is discussed below. ' 

The dependence of the observed profile shapes upon the magnetic 
field configuration and upon the properties of the atmosphere is so com- 
plicated that the only practical way of determining the fields is to cal- 
culate line profiles for a number of magnetic field strengths and 
inclinations and compare these to the observed profiles. The initial 
field strength derivation can be used to obtain starting values for the 
field strength, and profiles can be calculated for a given strength at a 
number of different orientations. Most theoretical procedures assume 
homogeneous fields. This assumption is a limitation in principle but 
not in practice, since the spatial resolution of the observations is 

srpall. The observed profiles represent averages over an area on the 

i 

Sun of approximately 175 by 1, 750 kilometers; they are averages over 
the depth range where the line is formed. The fields, therefore, could 
not be determined to smaller scale, even if the theory permitted it. 
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STARTING. POINT 


RESULT 



FIGURE 6. SUMMARY OF MAGNETIC FIELD CALCULATION 











The solution consists of intensities of the Stokes parameters; 
these can be combined to yield calculated intensities corresponding to the 
five observed polarization forms. The calculated profiles can then be 
compared to the observations, and the field strengths and inclinations 
in the sunspot will be assumed to be those which yield the best matches 
to the observed profiles. 

Because the theory is restricted to Zeeman triplets, this pro- 
cedure cannot be used to derive field properties from lines with more 
complicated Zeeman patterns. Another line in the observed spectral 
region has a large triplet splitting: the neutral chromium line at 
0. 5247574 micrometer (5, 247. 574 A), with Zeeman pattern ^ ^ . 

The same analysis could be performed for this line. Because the 
chromium line has approximately the same strength as the iron line 
(and, therefore, forms at approximately the same depths), its analysis 
would serve as a check on the iron line calculations but would not 
provide information concerning the fields at different depths in the 
sunspot. 
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INTRODUCTION 


The results of several studies of factors related to the practical 
use of solar magnetographs are collected in this section. Special 
emphasis is placed on uncertainties in the observed quantities and in 
their interpretation. The work reported was done between October 
1971 and June 1973. 

Several factors which influence magnetograph operation are 
discussed in this section, with special attention being given to uncer- 
tainties in the final results. A preliminary error estimate of ±20 
percent is suggested as a reasonable value on the basis of experience 
with other systems. 

The spectral line for which the Marshall Space Flight Center 
(MSFC) magnetograph was designed is shown to possess a number 
of disadvantages. Line profiles in sunspots cannot be calculated for 
this line, even approximately; the height of formation in the solar 
atmosphere is not known; and the strength of the line changes signifi- 
cantly for different regions on the solar disk. These disadvantages are 
largely a consequence of the fact that the line in question, A5250. 216 
of neutral iron, has a very low excitation potential (e. p. = 0. 121 
electron volt for lower level). A number of Zeeman triplet lines have 
lower excitation potentials of 2.0 electron volts or higher, and these 
lines would show these disadvantages to a much smaller degree. It 
is therefore recommended that serious consideration be given to 
choosing a new line of higher excitation as an alternative magnetograph 
line. Such a line also would be more likely to be observable at higher 
temperatures than those prevailing in sunspots and might permit the 
magnetograph to be used for study of magnetic stars. 



PROPERTIES OF THE NEUTRAL IRON LINE, X5250.216 

The MSFC Real Time Solar Magnetograph was designed to 
measure the degree of polarization in the neutral iron line at 525. 0216 
nanometers (5250.216 angstroms). This line has a lower-level excita- 
tion potential of 0. 121 electron volt; the central intensity in the photo- 
sphere as shown in the Utrecht Atlas is 0. 39 (Ref. 1); the equivalent 
width in Moore, Minnaert, and Houtgast's table of solar spectrum 
wavelengths is 62 milliang stroms (Ref. 2). The line is strengthened 
in stinspots. In a magnetic field, it shows a triplet Zeeman pattern 
with a strong splitting of about 4 milliangstroms per 100 gauss of field. 
The properties of the line are summarized in Table 1. 

The very low excitation potential of the lower level suggests 
that the line has a strong scattering component. That suggestion is 
supported by observations which show that the strength of the profile 
does not vary greatly with position on the solar disk; such behavior 
is characteristic of scattering lines. Therefore, a valid theoretical 
treatment must consider both the absorption and scattering mechanisms. 

Although a scattering matrix for the line has been derived 
(Ref. 3), the detailed calculation of a scattering line profile is very 
complicated. Obridko has done a calculation for a simple case (Ref. 4), 
but treatment of scattering alone using a realistic sunspot model is not 
yet practical. A calculation combining the effects of absorption and 
scattering for a realistic sunspot model is even more difficult. 

For this reason, most calculations of X5250 profiles assume 
pure absorption. Some recently published calculations of X5250 
profiles and of other profiles in sunspots are listed in Table 2. 
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TABLE 1. PROPERTIES OF Fe I A5250.216 


Transition 

a 5 D o - z 7 D, 

Lower and Upper 
Excitation Potentials (eV) 

0.121, 2.471 

Zeeman Pattern 

(0) 3 

Lands g Value 

3 

Log g/ 

(g is statistical weight of lower level; 

/ is oscillator strength for transition)* 

-4.46 

Photospheric Iron Abundance* 

log /NFe\ _ _ 5 2 

In I 

vw / 


*The proper choice of abundance and /-values of iron is not yet 
established; other published values are as justifiable as those 
given here. The gf-abundance product is better known than either 
of its constituent factors. 
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TABLE 2. LINE PROFILE CALCULATIONS IN SUNSPOTS 


AUTHOR 

SPECTRAL LINE 

SUNSPOT MODEL 

PHOTOSPHERIC MODEL - 

REFERENCE- 

Moe and Maltby 

X5250.2 


Holweger 

5 

Hlnoux 

Na D, Mg b (wings) 

Henoux 

- 

6 

H§noux 

Weak to medium Fe I, Fe II, 
Cr I, Cr II, Ti I, 

Ti II 

Henoux 


6 

Yun 

Na D (wings) 

Yun 


7 

Beckers 

X5250.2 

■ 


Bilderberg Continuum 
Atmosphere- 

8 

Gohring 

X5250.2 


Arbitrary n 0 

9 

Gfohring 

X5250.2 


Holweger 

9 

Moe 

X5250.2 


Combination of 
Muller-Mutschlecner 
and Goldberg 

10 

Staude 

X5250.2 

H§noux 
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Although some of this work was elaborate, only one author, Staude, 
has published profiles of A.5250. 2 calculated with a sunspot model. 

(Moe states that he has succeeded in making a similar calculation 
but has published no results. ) As will be shown later, the profiles 
calculated in sunspots do not appear to be valid. 

In sunspot spectra, a weak line of titanium oxide is blended with 
the red wing of A5250, and this blend should be considered in a precise 
profile calculation. Each line modifies the radiation field seen by the 
other, and a rigorous approach requires that the source functions for 
both lmes be modified correspondingly. Since in this case the strengths 
of the lines are so different, it is probably acceptable to assume 
that the titanium oxide line has negligible effect on the source function 
for the iron line and that the titanium oxide line is so weak that the 
effect of the modification of its source function by the iron line would 
be found to be of the order of the observational uncertainty or smaller. 
Then the separate profiles can simply be superimposed; the residual 
intensity for the combined line will be given by: 

' S(X >Fe+TIO = R(X) Fe ' t 1 ' 00 ' R <*W' 

No blend calculations, even on this level of approximation, have been 
published. Moe has calculated but not published sunspot profiles for 
X5250 in which he treated the presence of blends by adjusting the con- 
tinuum level. His reasoning is that the apparent continuum under the 
conditions prevailing in sunspots is likely to be crowded with weak, 
unresolved molecular lines. This approach yields profiles which 
agree well with observations, but the physical validity of the argument 
is questionable. Moreover, it will be shown that in this case the fact 
that calculated profiles agree with observations is not sufficient 
evidence that the calculation is correct. 
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The low excitation potential of the lower level of the line causes 
the absorption coefficient to be sensitive to the temperature and electron 
pressure in the model atmosphere. Figure 1 shows the fraction of 
iron atoms in the lower level as a function of depth in the sunspot model 
atmosphere (the model is that of H^noux, Ref. 6). The top curve in 
the figure, marked by triangles, is the relation for X5250. The lower 
curve, shown for reference, is the corresponding set of values for the 
neutral iron line at 617.3348 nanometers. This line has a lower 
excitation potential of 2. 213 electron volts (Ref. 1.2); it is used for 
comparison because self-consistent profiles can be calculated for it 
in sunspots. The number of atoms capable of absorbing X6173 decreases 
by about a factor of ten with increasing height in the atmosphere, but 
the number for X5250 is nearly constant over most of the range and 
is close to its peak value at the top of the atmosphere. [ The values 
were calculated using the Boltzmann and Saha equations, assuming that 
iron could exist in neutral and singly-ionized forms; this calculation, 
of course, implicitly assumes the atmosphere is in local thermodynamic 
equilibrium (LTE). ] 

Figure 2 shows q , the ratio of line -to -continuous absorption 

A 

coefficients, as a function of depth in the H^noux model atmosphere. 

Again, the triangles are used to signify X5250 and circles to signify 

X6173. The values for both lines rise sharply toward the top of 

the atmospheric model, primarily because the opacities fall rapidly 

over the same range. The effect of the increase in r\ is much more 

A 

serious for X5250 because of the; different line absorption coefficients 
for the two lines; the coefficient for X6173 is much smaller at the top 
of the atmosphere, compared both to its own peak value and to the 
values for X5250. The coefficient for X5250 declines more slowly 
with height and has a much larger average value. Figure 2 also 
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NOTE: 

EQUALS THE FRACTION OF 
IRON ATOMS IN THE LOWER 
LEVEL OF THE LINE-FORMING 
TRANSITION. * - A5250.216, 
o = A61 73.348. 




log n 




provides the reason why Moe's addition of continuous absorption 
improved the appearance of his calculated profiles without correcting 
the fundamental difficulty: adding continuous absorption will reduce 
the ordinate values of the curves of Figure 2, but the values of the 
line absorption coefficients also are important to the result, and they 
are unchanged by the addition of continuous absorption. 

When combined, these effects are sufficient to make impossible 
a self-consistent LTE profile calculation of X5250 in sunspots. The 
residual intensities, R. , may be evaluated by direct integration of 






ratio of line -to -continuous absorption at 
wavelength X . 
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The calculations were done at the disk center (p = 1.00) for 
zero magnetic field. (The addition of a magnetic field will not affect 
the derived contribution relation, since the line -to -continuous absorp- 
tion ratio for circularly polarized light in a longitudinal magnetic field 
is equal to the value in the zero field case. ) The ^-values were taken 
from Corliss and Warner (Ref. 13). The iron abundance (log j© Np^/N^j) 
was set equal to -5. 20 (Ref. 14). The microturbulent velocity in the 
sunspot was assumed to be constant with depth and equal to 1. 0 kilo- 
meters per second. 

The contributions of the various depths to the emergent intensity 
are displayed in Figure 3 by plotting the differences, AR^/A (log 1Q T 0 )» 
as a function of log r ; the A R> values were obtained from the 
integration of Equation (2). The A (log j 0 t 0 ) intervals in each case are 
the depth intervals at which the H£noux model is tabulated; the differences 
are plotted at the centers of their depth intervals. 

As Figure 3 shows, the intensity contribution to X5250 is at its 
maximum at the top of the model, with the peak value being undefined. 
Integration of these values over depth clearly will not yield a completely 
defined emergent intensity. On the other hand, the contributions to 
X6173 are essentially all within the depth range of the model, so the 
residual intensity obtained by integrating those values over depth can 
be considered to be significant. 

For photospheric conditions, the temperatures and electron 
pressures are high enough that the number of atoms in the lower level 
of the X5250 transition is reduced, and self-consistent LTE profiles 
can be calculated. The X5250 profiles calculated for photospheric 
models are therefore valid insofar as their underlying assumptions 
are valid. However, it does not seem reasonable to expect that profiles 
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calculated for photospheric models will also be representative of sunspot 
conditions, and they therefore should not be used to interpret sunspot 
observations. (In practice, most magnetographs which use the X5250 
line and depend on calibration from calculated profiles were calibrated 
using profiles calculated with photospheric models. ) 

Because X5250 is quite strong in sunspots, its profile is often 
saturated. When saturated profiles in different polarization states 
are superimposed, the intensity distribution of the combination may 
appear quite different from the distributions in the separate components. 
Often, in fact, the combination appears to have a spurious n component, 
usually shifted somewhat from the normal position (Refs. 5, 9, 11, 15, 
and 16). This so-called "anomalous n component" further complicates 
the analysis of X5250 observations. 

Another source of complication of the theoretical treatment of 
X5250 is the effect of polarization of the atomic sublevels involved in 
the transition. If the splitting in a magnetic field is not complete, so 
that different sublevels overlap, then the sublevels are not independent, 
and phase relations exist between their wave functions. This inter- 
ference between sublevels is called "level-crossing interference". 

The Hanle effect is an example of level-crossing interference. 

The sublevel polarization is distinct from the polarization 
introduced by the Zeeman splitting. unusually brief and clear 
description of the various sources of polarization has been given by 
Lamb (Ref. 17); it is quoted in the following paragraph. 

The polarization of light in solar absorption and 
emission lines may result from the action of one or 
more of three distinct processes. First, in the pres- 
ence of a magnetic field the light in absorption or 


54 



emission lines may become partially polarized simply 
as the result of Zeeman splitting. Although the radia- 
tion absorbed or emitted at each frequency in the line 
by an assembly of atoms is then at least partially 
polarized, if atomic level polarization does not accom- 
pany the Zeeman splitting and if the assembly is optically 
thin in the line the total radiation absorbed or emitted, 
when integrated in frequency over the whole of the line, 
will be isotropic and unpolarized (one speaks of no "net" 
polarization of the line). Second, the light in Zeeman 
split absorption or emission lines formed in an assem- 
bly which is not optically thin may show a different 
partial polarization at each frequency from that of an 
optically thin assembly. Even in the absence of atomic 
level polarization, this phenomenon usually leads to net 
polarization of the light in the line. Finally, light in 
absorption or emission lines may become polarized as 
a result of the polarization of one or both of the atomic 
levels involved in the formation of the line, since in this 
case the assembly will preferentially absorb, emit, and 
scatter radiation of a particular polarization and angular 
distribution in radiative processes beginning at each of 
the polarized levels. Generally speaking, atomic level 
polarization leads to net polarization of the light in the 
line, whether or not the line also undergoes Zeeman 
splitting or the assembly is optically thin. 


The study of these effects is still in its infancy; both Lamb and 
L. L. House have been working actively in this field (see Ref. 18, for 
example). A fairly complete analysis can be done only for certain 
simple resonance lines. The sublevel structure of A.5250 is far too 
complicated to permit an analysis of that transition with the techniques 
now in use. It should be observed, however, that Lamb (Ref. 17, 
p. 159) has estimated from collisional relaxation rates that "at photo- 
spheric densities the energy 'overlap' of states within the same atomic 
level will not be significant [for A.5250] when the magnetic field strength 



is greater than 10 2 G. M Presumably, level-crossing interference also 
will be unimportant for lines formed in the lower depths of sunspots. 

For lines formed at greater heights, however, level-crossing inter- 
ference may become important. The exact height of formation of A.5250 
is uncertain, but it probably is formed in the upper photosphere. There- 
fore, the possibility of the existence of level-crossing interference in this 
line cannot yet be eliminated. 



ALTERNATIVE SPECTRAL LINES FOR 
MAGNETOGRAPH OPERATION 

Because a theoretical calibration from calculated sunspot pro- 
files for \5250 cannot now be accomplished, it may be advantageous to 
consider adopting another spectral line for magnetograph operation. 
With this in mind, a number of Zeeman triplet lines within about 200 
angstroms on either side of X5250 were examined for suitability as 
magnetograph lines. The results of this preliminary study are sum- 
marized in Table 3. The lines were first checked for the presence of 
blends, using the photospheric tracings in the Utrecht Atlas (Ref. 1) 
and a photographic sunspot spectrum borrowed from the Kitt Peak 
National Observatory. The lines which did not appear to be badly 
blended in those two references were investigated further. In Table 3, 
M UA" signifies Utrecht Atlas; "KPNO M , the Kitt Peak photographic 
spectra; and "MMH", the Moore, Minnaert, and Houtgast solar wave- 
length table (Ref. 2). The spectral line identifications in this last 
reference are useful for detecting and identifying blends. 

It must be emphasized that probably no strong line can be found 
to be entirely free of blends in sunspots because of the large number of 
molecular lines which can form under these conditions. The presence 
of weak blended lines therefore cannot disqualify a spectral line for 
magnetograph use. 

Spectral lines at greater distances from X5250 also might be 
considered. Some interest now exists in the polarization of lines in 
the near infrared (Ref. 19)> but these would not be suitable for use in 
a magnetograph located near sea level because of the strong telluric 
absorption by water vapor at long wavelengths. Harvey (Ref. 20) has 
published a list of lines with large Zeeman splittings, which could be 
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TABLE 3. PRELIMINARY SURVEY OF SPECTRAL LINES NEAR 5200 ANGSTROMS 


WAVELENGTH (A) 

RESULT OF BLEND STUDY 

CONCLUSION 

5082.351 

KPNO and UA spectra suggest blends are not serious; according to MMH, 

Probably best avoided because of 

(Nil, multi pi et number 130) 

it may be blended in the wings with lines of C 2 and MgH 

• gX 2 = 39 x 10 _1 ° cm 2 (von Kluber) 

• Transition: z *P° - e *P t 

• Excitation potential of lower level = 3.642 electron volts 
e According to von Kluber, line is weakened in sunspots 

small splitting and weakness in 
sunspots 

5145.104 

Both KPNO and UA spectra show serious blends 

Not usable 

5202.351 

KPNO spectrum shows serious blend, confirmed by MMH 

Not usable 

5215.190 

UA spectrum shows blend in red wing; continuum is depressed in this 

Not usable 

(Fel, multi plet number 553) 

region. 

• Transition: z S D° - e s Dj 

• Excitation potential of lower level = 3.26 electron volts 


5217.398 

Profile looks fairly clean in UA and KPNO spectra; MMH suggests 

Probably could be used but is not 

(Fel, multiplet number 553) 

■ ' ' ' ■ ' ■ ■ ' ' 

possible blend 

• gX 2 = 41 x 10‘ 10 cm 2 

• Transition: z S D° - e s D t 

• Excitation potential of lower level = 3.197 electron volts 

a good choice because of small 
splitting (see entry for X5263) 

5229.862 

Profile looks clean in UA but KPNO spectrum shows possible blend 

Probably could be used but is not 

(Fel, multiplet numbers 553, 1090) 

in red wing 

a good choice because of small 

• gX 2 =41 x 10" 10 cm 2 

• Transition: z s Dj - e 5 D« 

• Excitation potential of lower level ■= 3.269 electron volts 

splitting (see entry for X5263) 

5247.576 

• gX 2 = 69 x 10" JD cm 2 

This line has no advantages over 

(CrI, multiplet number 18) 

• Transition: a 5 D 0 - z 5 P? 

X5250 


• Excitation potential of lower level = 0.957 electron volts 

This line, like X5250, would be an excellent magnetograph line if it 
had a higher excitation potential. Attempts to calculate absorption 
profiles for sunspot models have shown that a profile cannot be 
defined for the same reasons that profiles are not defined for X5250. 





TABLE 3 


Concl uded 


WAVELENGTH {A°) 

RESULT OF BLEND STUDY 

CONCLUSION 

5253.470 

(Fel, multi plet number 553) 

Appearance of profile is good in UA and KPNO spectra; MMH shows blends 
with two faint lines 

• gX 2 = 41 x IQ* 10 cm 2 

• Transition: z S D? - e s Di 

• Excitation potential of lower level = 3.27 electron volts 

Splitting may be too small for 
line to be a good choice. It is 
probably usable. 


Zeeman splitting is not detectable in the KPNO spectrum. For a field 
of 2000 gauss, the value of AX^ should be about 38 milliangstroms, 
compared to about 90 milliangstrom for X6173 for about the same field 
strength. If a AX^ as small as 2.5 milliangstroms could be resolved, 
the lower limit on measurable field strength would be about 130 gauss. 


5263.316 

Appears badly blended in UA and KPNO spectra 

Not usable 

5273.172 

Appears badly blended in UA and KPNO spectra 

Not usable 

5283.631 

Appears badly blended in UA and KPNO spectra 

Not usable 

5302.309 

(Fel, multi plet number 553) 

Profile looks clean in UA but KPNO spectrum appears blended 

• gX 2 = 42 x lo- 1D cm 2 

• Transition: z 5 D° - e s Dz 

• Excitation potential of lower level = 3.269 electron volts 

Not usable because of conspicuous 
blend in sunspots and small 
splitting 

5324.193 

(Fel, multiplet number 553) 

Profile looks clean in UA but appears to be blended with a weak line, 
possibly Cr I, in KPNO 

• gX 2 = 42.5 x IQ" 10 cm 2 

• Transition: z 5 D° - e s Di, 

• Excitation potential of lower level = 3.20 electron volts 
For a field of 2000 gauss, AX^ = 39.6 milliangstroms 

Probably usable. {This is the line 
used by the Aerospace Corporation 
magnetograph.) 

5339.939 

Profile looks blended in both UA and KPNO spectra 

Not usable 

5393.178 

(Fel, multiplet number 553) 

Profile looks blended in both UA and KPNO spectra, confirmed by MMH 

Not usable 




examined. Zwaan and Buurman have compiled a list of sunspot lines 
which are very weak in the photosphere and which therefore should not 
be distorted in shape by the effects of stray light. This list does not 
appear to have been published. One of the lines listed, according to 
Reference 21, is X6064. 626 of neutral titanium. Another relevant 
reference is Reference 22, in which Wittmann reports the results of 
examining a number of red lines for the presence of blends. This 
material would form a starting point for a more extensive search of the 
solar spectrum. 



SOURCES OF ERROR FOR SOLAR MAGNETOGRAPH 


EFFECTS WHICH DEPEND SOLELY ON POSITION ON THE 
SOLAR DISK 

As is well known, the radiation observed at the center of the 
solar disk characterizes a greater range in depth than that observed 
at the limb. A crude numerical estimate of this effect can be made by 
using the Eddington- Barbier relation, 

Tq = cos 9 = (i , 

where Tq represents the maximum optical depth in the solar atmosphere 
which can be observed at solar position angle 0 . 

Since the geometrical depth, d, in the atmosphere is approxi- 
mately proportional to the log of the optical depth, 

d ph - K p h log T ph 

d pe “ K pe lQ g T pe 
d u - K u log T U r 

where 

ph - photosphere 
pe - penumbra 
u - umbra. 


61 




In general, 


Kph ^ ^ ^u 


and 


^ph ^ ^P® ^ ^ * 


The approximate form of the relation between d and p across the 
solar disk can be shown by letting K = 1. The zero point in the geo- 
metrical depth scale is arbitrary; here it has been chosen so that 
d = 1 at p = 1, which is the deepest level observed. 

d 

1,0° n 
0.955 
0.903 
0.846 

0.779 Increasing 

0. 700 depth 

0. 603 
0.478 
0. 301 
0.000 

Thus, for a given solar feature the depth observed at p = 0. 5 
will be about 0. 7 of that observed at p = 1. This property can be very, 
useful, since it allows a partial determination of the variation of observ- 
able quantities with depth in the Sun. 

Foreshortening toward the limb will tend to decrease the spatial 
resolution of the observations. The magnetograph aperture has a pro- 
jected size of about 5 minutes of arc square, at the distance of the Sun 
from the Earth. The aperture will accept a distance on the Sun at the 


I 

1.00 
0.90 
0.80 
0.70 
0.60 
0.50 
0.40 
0.30 
0 . 20 
0. 10 


62 



center of the disk of 2. 17 X 10 s kilometers. If the aperture is placed 
so that its outer edge is just at the limb, the linear distance it subtends 
on the Sun becomes 5. 84 X 10 5 kilometers, or more than two and one- 
half times as large as at the disk center. 

Finally, any observation using line profiles will be affected to 
some extent by the tendency of profiles to broaden and flatten toward 
the solar limb. This tendency should not be reflected in the degree of 
polarization in the profile, so the effect for a polarization magnetograph 
will be mostly that of reducing the signal level somewhat, since both 
the brightness of the Sun and the depth of the line profile, which provides 
contrast, are reduced. The "limb darkening" is observed to be smaller 
in the umbrae of sunspots than in the photosphere; the photospheric and 
penumbral limb darkening are approximately the same. 

EFFECTS CAUSED BY CHANGES IN STRENGTH AND SHAPE OF 
LINE PROFILE 

One assumption which is implicit in the entire magnetograph 
design and operation is that the line absorption coefficient in each of 
the Stokes parameters is identical to the line absorption coefficient in 
the absence of a magnetic field, although it may be shifted in wavelength 
by an amount which corresponds to the Zeeman splitting. Knowledge of 
the form of the line absorption coefficient is, in principle, sufficient to 
determine the shape and strength of the line profile corresponding to a 
transition between known levels in a model atmosphere of known proper- 
ties; therefore, the interpretation of the magnetograph results can be 
derived from the properties of the line profile thus determined. It 
follows that anything which alters the shape cr strength of the line pro- 
file on the Sun will affect the interpretation of the magnetograph, 
measurements. 
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The profile of \5250.216 is particularly susceptible to altera- 
tion by varying local conditions on the Sun. Since determining and using 
a separate magnetograph calibration for each of these local areas is 
difficult, an average calibration usually is used. The local variations 
will then result in errors when the average calibration is used to inter- 
pret results in the variant regions. 

The local variations can be dependent on time and on position 
of the solar disk. They also vary in magnitude. A quantitative estimate 
of their effects is almost impossible to achieve, but qualitative descrip- 
tions can be given. 

The effect on line profiles of motion of the solar gas is well 
known. If the scale of the motion is small compared to the range of 
depths over which the lines are formed (microturbulence), the line 
profile is broadened and its equivalent width is increased. If the scale 
of the motion is large compared to the range of line- formation depths 
(macroturbulence), the profile is broadened but its equivalent width 
is unchanged. Large-scale motions may also shift the entire profile 
in wavelength or cause it to become asymmetrical. The turbulence 
effects depend upon position on the solar disk: the profiles tend to 
broaden and flatten toward the limb. 

The normal variation toward the limb can be determined and 
compensated for in the reduction and analysis. (If it is not compen- 
sated, errors of residual intensity of 0. 10 or more -- 10 percent of 
the continuum value -- can result. ) In a system which uses a filter to 
isolate the line, and does not have an entrance slit, shifts of the entire 
profile must be compensated in that way also. Most of the uncompensable 
error will come from undetected asymmetries and from unresolved 
random motions in umbrae, penumbrae, and active regions. The 
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uncompensable errors in residual intensity in each Stokes parameter 
probably will be less than 0.05; the relative error in each measurement, 
of course, will depend upon the part of the profile in which the measure- 
ment is made. Since the average degree of circular polarization, 
is found by combining signals measured in two different polarization 
states, and since both signals may contain errors which are comparable 
in absolute value and partially interdependent, the resulting error in 
P v cannot be estimated for a general case. Wiehr (Ref. 23) discusses 
errors caused by mis centering of the line profile in the original polarim- 
eter designs used at the Crimea observatory and at the Locarno station 
of the Gottingen observatory; these errors can be very large. 

Of the spectral lines used for magnetic field determinations, 
\5250. 2 shows some of the largest natural variations in strength. The 
line has moderate strength in the photosphere; it is stronger in sunspots 
and weaker in photospheric magnetic regions. Chapman and Sheeley 
(Ref. 24) give the photospheric equivalent width as 62 miiliang stroms; 
the equivalent width they measured in photospheric magnetic regions is 
27 percent smaller. The equivalent width in sunspots is difficult to 
determine because the spot spectra are usually contaminated by stray 
light, and no reliable estimate has been found. 

The strength variations are thought to be caused partly by 
temperature effects and partly by magnetic saturation (Refs. 24, 25, 
and 26); the temperature effects appear to be confined to the line cores. 
Whatever the cause, the variations in profile strength have been 
responsible for discrepancies between magnetograph results obtained 
with X5250 and with other lines. Usually, the calibration for X.5250 
is based on an "average" profile for the solar disk and thus cannot 
account for the strength variations. Harvey and Livingston (Ref. 27) 
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found that the fields they measured using X5250 appeared to underesti- 
mate the true field strength by a factor of about 2 in the linear portion 
of the profile and about 5 near the line core; Chapman and Sheeley 
(Ref. 24) have estimated that a calibration using an average disk profile 
of X5250 will yield magnetic field strengths which are too low by a factor 
of 1.5 to 2; and measurements made with molecular lines also differ 
consistently from the X5250 values (Ref. 28). 

Therefore, calibrations based only on disk-average profiles 
of \5250 generally appear to yield estimates of magnetic field strength 
which are much too low. Some supplementary information is needed. 

Indeed, Harvey et al (Ref. 25) now use data from X5233 to establish the 
average field strength in the areas they observe, and use the X5250 
observations to provide finer detail. 

ERRORS WHICH ORIGINATE WITHIN THE TELESCOPE AND 
THE MAGNETOGRAPH 

The magnitudes of errors arising within the magnetograph and 
the telescope will be determined by making measurements of auxiliary 
sources or solar regions whose characteristics are known. Although 
a complete description of the errors cannot be made until that has been 
done, a general review of some of the possible sources of error can 
be given here. 

The magnetograph contains a flat glass plate to be used to com- 
pensate for residual polarization introduced by the telescope. Even after 
this compensation has been made, a careful check should be made for 
residual polarization in all modes to be measured, because residual 
polarization is very difficult to remove completely. Furthermore, if 
any residual polarization can be measured, it may be possible to remove 
its effect in the data analysis. Residual polarization studies should be 
made for a number of different source intensities and telescope orientations 
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One of the chronic difficulties experienced with umbral and pen- 
umbral observations is contamination of the fainter solar intensities 
by stray light from brighter regions. The scattering processes which 
redistribute the light occur in the Earth’s atmosphere and within the 
telescope. The net polarization produced by the atmospheric scattering 
should be negligible, but the redistribution of light at the telescope 
optical surfaces may introduce a net polarization. Thus, a small 
polarized component of stray light may exist which would tend to dis- 
tort the relative intensities measured in different polarization states. 

This component, if present, would probably be essentially linear and 
therefore less serious for the circular polarization measurements. 
Nevertheless, it is a possibility that should be examined, since such 
polarized stray light might cause an appreciable distortion if the signal 
level is very low. 

A more serious possibility is that of intercommunication or 
’’crosstalk" between signals. If, for example, some fraction of unpolar- 
ized intensity is mixed with a circularly polarized component, a sizable 
error can result. Instruments which use electrically switched polarizing 
crystals are subject to this type of signal leakage, so the possibility of 
its existence should be examined (see, for example, Reference 25). 

ERRORS CAUSED BY LIMITED RESOLUTION 

The Cassegrain telescope to which the magnetograph is 
attached has a diameter of about 30, 5 centimeters. The corresponding 
Rayleigh criterion value at \= 5250 angstroms is 0. 43 second of arc; this 
will be assumed to be the diffraction-limited spatial resolution of the tele 
scope. The design goal for the magnetograph optics is a spatial resolu- 
tion of 2. 5 seconds of arc for the 5-by-5 minute of arc field size, and 
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1 second of arc for the 2-by-2 minute field size obtained with the 
magnification system (Ref. 29, p. 2-11). The resolution values 1 or 
2,5 seconds thus represent the best possible performance of the system; 
the resolution may be further limited by atmospheric seeing. 

The typical daylight seeing values at the telescope tower have 
not yet been completely determined; in most locations, however, 1 
second of arc is considered excellent, and 2 seconds or more are more 
usual (Ref. 12, p. 19). Moreover, the seeing usually varies rapidly 
with time, even on the best days, so an observation of more than a few 
seconds in length probably will show distortion. A reasonable predic- 
tion would be that on average days, observations made with the 5-minute 
field size should not be severely seeing -limited; the full resolution of 
the 2-minute field will be attained only at times of exceptionally good 
seeing. 

Good spatial resolution is required for precise mapping of 
magnetic fields on the Sun. (A really fine discrimination of position 
is not possible from the Earth, since here an angle of 1 second of arc 
corresponds to a linear distance of about 725 kilometers at the center 
of the solar disk. ) Spatial resolution is also important because the Sun 
shows many surface inhomogeneities of sizes near the limit of res- 
olution or smaller. At maximum resolution, the smearing of these 
inhothogeneities will be smallest, and the resulting distortions of the 
observations will be minimized. Several of these inhomogeneities are 
briefly described below: 

• Umbral dots are small, bright transient regions 
which appear in sunspot umbrae. The mean diam- 
eter of the dots measured by Beckers and Schroter 
(Ref. 30) was 420 kilometers (about 0. 6 second of 
arc); the mean lifetime was 1500 seconds. (The life- 
time was defined as the time between the occurrence 
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of half-maximum brightness in the increase and 
decrease of the dot's intensity. ) Umbral dots may 
be associated with inhomogeneities in the umbral 
magnetic field. 

• Penumbral filaments (Ref. 12, p. 75) are irregular, 
elongated bright regions in the penumbra. The lengths 
of the filaments vary a good deal; a representative 
value for large spots is about 10 seconds of arc. Their 
lifetimes are of the order of a few hours, so a promi- 
nent penumbral filament in an observed region probably 
would be noticed. 

• Magnetic knots are small regions in the photosphere 
which are characterized by unusually strong magnetic 
fields (Ref. 31). The strengths of many spectral lines, 
including \52 50, decrease in these regions. Atypical 
field strength is 1000 gauss; diameters are about 1000 
kilometers. Lifetimes of knots are greater than 30 
minutes. 

• Pores are small sunspots with no penumbrae (Ref. 31). 

Their diameters are 2 to 5 seconds of arc, so at times, 
of very poor seeing they may not be resolved. They 
typically have brightnesses of about half that of the sur- 
rounding photosphere and field strengths greater than 
1500 gauss. 

• Nonmagnetic gaps (Ref. 31) are regions in the photo- 
sphere where the strength of certain spectral lines 
decreases and the continuum brightness increases, 
but there is no strong local magnetic field. The diam- 
eters of these regions are about 1100 kilometers, and 
lifetimes are about 40 minutes. 

Any of these inhomogeneities is capable of producing some 
distortion of the magnetograph record, especially at times of poor 
seeing. The amounts of such distortion will be determined largely by 
the magnitudes of the effects of errors arising within the telescope and 
magnetograph. 



It is apposite to quote here a statement made by Deubner and 
Liedler in their discussion of the magnetograph at the Fraunhofer 
Xjjistitut (Ref. 32): "In attempting to calibrate the vectormagnetograph, 
it has been shown here that due to our limited knowledge of the nature 
of these inhomogeneities (on the solar surface), which doubtless exists, 
the construction of a universal calibration curve is more or less impossi- 
ble or even meaningless. Such a calibration would be valid only for the 
single spot from which it was derived during a particular state of develop- 
ment, and only for the prevailing spot field. " 

ERRORS ASSOCIATED WITH UNCERTAINTIES IN THEORY 

Calibration curves relating magnetograph signal levels to 
magnetic field strengths on the Sun must be based, ultimately, on 
some theory of line formation in magnetic fields. As will be shown 
later, theory does not yet permit a precise prediction of signal levels 
for all conditions on the Sun. 

Approximate calibration curves have already been calculated. 
While they are not exactly applicable, they should serve to show the 
general behavior of signal as a function of field strength. The rela- 
tions thus obtained consist of a series of nested curves, one for each 
value of the field orientation angle. These curves are fairly closely 
spaced; since in some cases the spacing may be of the order of the 
uncertainty of the signal, the interpretation could be ambiguous even 
supposing that the calibration curves were precisely known. More- 
over, the curves are not monotonic; the signal increases with increas- 
ing field strength until a field of the order of 1000 gauss is reached; as 
the field increases beyond that value, the signal drops abruptly. There- 
fore, an additional ambiguity exists in the region around the peak unless 
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independent information is available to show on which branch of the 
curve the signal lies. To make matters worse, at present there is no 
way to determine an absolute calibration curve for \5250, so given signal 
values cannot be attributed to specific field strengths. This comes 
about partly through the inherent difficulties of calculating the polarized 
intensities in a solar magnetic field, and partly through the properties 
of the transition producing X5250.216. 

Only one or two transitions can be treated in detail in solar 
magnetic fields, and the calculations for these cases are so cumber- 
some as to be impractical for the large number of repetitions needed 
to define a calibration curve. House and Cohen (Ref. 33) described the 
calculations necessary to define the scattering in a magnetic field (with 
simplifying approximations) for a resonance line which gives rise in a 
magnetic field to a normal Zeeman triplet. The field must be weak. 

The transition considered is 1 S 0 -♦ 1 P 1 , which for the solar astro- 
physical case would be found only in the resonance lines of neutral 
calcium and magnesium, in the visible wavelengths. 

The 5250. 216 angstrom line comes from the transition 5 Do *" 7 Dj 
(R ef<i 34); the lower excitation potential is 0. 12 electron volt, the upper 
2.47 electron volt. The line is therefore one of low excitation but not a 
resonance line. One would expect it to have a strong scattering component, 
and this expectation is reinforced by its behavior on the Sun, where the 
variation in line strength from center to limb is less than would be 
expecred for a line formed in pure absorption (Ref. 4). Nevertheless, 
the line also has a strong absorption component. It cannot be correctly 
treated by either a pure scattering or a pure absorption approximation, 
and these are the only computational methods which are practical for the 
large number of calculations required. The two most often used pro- 
cedures for calculations in a magnetic field, described in References 10 
and 35, require that pure absorption or pure scattering be assumed. 
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There is general agreement that, because of the intrinsic differ- 
ences between the regions, calibration curves established for the 
photosphere should not be used in interpreting sunspot observations. 

The calculation of expected signal values in sunspots is even more 
difficult than in the photosphere. All the difficulties of the photospheric 
calculation remain, and in addition the existing sunspot models do not 
permit a self-consistent profile calculation under any set of assump- 
tions (Ref. 36). 

The interpretation of the observations is even more uncertain 
because of several magneto -optic effects which can alter the plane 
of polarization after it is produced in the spectral line transition on 
the Sun. The Hanle effect usually is mentioned in this connection; 
others may enter as well (see Ref. 37, for example). The effects are 
much more severe for linear than for circular polarization. Indeed, 
Beckers has stated, "There is no simple relationship between the 
amount and direction of linear polarization and the strength and direc- 
tion of the magnetic field. Only with very large approximation can one 
say that the degree of polarization is related to the total field strength 
and the direction to the longitudinal magnetic field strength. " (See Ref. 
38, p. 13.) 


Many magnetograph calibrations take advantage of the fact that, 
for small fields, the measured signal is (approximately) directly pro- 
portional to the field strength. This relation can be used if it is under- 
stood that it applies to small fields only and that the proportionality 
"constant" depends explicitly on position on the solar disk, so that the 
value used must be correct for the solar coordinates of the sunspot 
(see Ref. 12, p. 194, and Ref. 29, p. 1 — 17). Adjustment of the 
"constant" according to disk position will be discussed in a later sec- 
tion of this report. 
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PROBABLE ERRORS IN FIELD STRENGTHS DETERMINED FROM 
ZEEMAN SPLITTING 

The results of a determination of magnetic fields in sunspots 
from the Zeeman splitting of X 5250 illustrates the error which can 
enter a well-documented field strength determination/ The observations 
and method of deriving field strengths are described in Section 1. 

The calculated field strengths were reviewed to estimate the error in 
the final values. Most of the field strengths are averages of several 
values obtained from separate spectra made of the same location in 
the sunspot, so the deviations of the separate values from the average 
form a part of the error estimate. When the calculated field strength 
values were reviewed, the identification of the separate Zeeman com- 
ponents was also reviewed. In some cases, the identifications were 
revised, causing changes in the average field value for the Location. 

These changes generally were small. A value which departed greatly 
from the average of the other values was discarded in one location 
(square Xl-g of the magnetic field map). The final map is shown in 
Figure 4; it does not differ appreciably from the one given in Reference 
39. The deviations from the mean values shown in the map are tabulated 
for each location in the sunspot in Table 4. 

The average deviation for all the values from their separate 
means is 180 gauss; the uncertainty caused by scatter in the individual 
magnetic field values for this sunspot, therefore, is ±180 gauss. 

The final values are uncertain for other reasons also. Because 
the splittings are determined from the observed positions of the Zeeman 
components, uncertainties in wavelength will affect the results. The 
spectral intensities were recorded at wavelength intervals of 5. 88 milli- 
angstroms, so this interval size imposes a limit to the wavelength resolu- 
tion. If the location of a Zeeman component is defined as the wavelength 
of its minimum measured intensity, the uncertainty in that location is at 
least ± 5. 88/2 milliangstroms. 
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FIGURE 4. REVISED LINE-OF-SIGHT FIELD STRENGTHS IN GAUSS 
FOR SUNSPOT OBSERVATIONS OF JULY 14, 1969 











TABLE 4. DEVIATIONS OF SEPARATE FIELD STRENGTH VALUES 
FROM MEAN VALUE FOR LOCATION 


LOCATION IN 
SUNSPOT 

MEAN FIELD 
STRENGTH (gauss) 

NUMBER OF 
VALUES AVERAGED 

AVERAGE DEVIATION 
FROM MEAN (gauss) 

Ilia 

908 

2 

152 

mb 

805 

3 

67 

IIIc 

605 

1 

0 

illg 

906 

1 

0 

IVa 

1,250 

4 

609 

IVb 

1,230 

3 

686 

IVg 

1 ,070 

5 

315 

IVh 

830 

2 

75 

I Vi 

151 

1 

0 

V, Via 

1 ,680 

7 

155 

V, VI b 

1 ,670 

7 

139 

V, Vic 

1,420 

9 

279 

V, VId 

604 

2 

151 

V, VIg 

1,430 

10 

300 

V, VIh 

1,330 

3 

419 

V, VI i 

944 

2 

416 

VII, Villa 

1,780 

5 

126 

VII, VUIb 

1,670 

8 

140 

VII, VI I I c 

1 ,360 

4 

150 

VII, VII Ig 

1,900 

6 

120 

VII, VHIh 

1,400 

4 

376 

VII, VII Ii 

926 

4 

284 

IX, Xa 

2*060 

5 

104 

IX, Xb 

1 i770 

5 

313 

IX, Xc 

1*640 

10 

211 

IX, Xd 

1 ,360 

1 

0 

IX, Xg 

1*990 

5 

87 

IX, Xh 

2,010 

5 

206 

IX, Xi 

1,680 

5 

268 


75 






TABLE 4 


Concluded 


LOCATION IN 
SUNSPOT 

MEAN FIELD 
STRENGTH (gauss) 

NUMBER OF 
VALUES AVERAGED 

AVERAGE DEVIATION 
FROM MEAN (gauss) 

XIa 

1,940 

5 

126 

Xlb 

1,680 

5 

139 

XIc 

1,300 

4 

169 

Xlg 

2,060 

4 

58 

Xlh 

2,020 

5 

202 

XI i 

1,510 

5 

120 

Xlla 

1,750 

5 

277 

XI lb 

1 ,450 

5 

313 

XI Ic 

1 ,440 

2 

225 

XI Ig 

1 ,840 

5 

134 

XI Ih 

1,670 

5 

131 

XI I i 

1,120 

4 

170 

XI 1 1 a 

1,210 

5 

392 

X 1 1 1 b 

755 j 

2 

0 

XI I Ig 

1 ,580 

5 

360 

XII Ih 

1 ,240 

5 

187 

XI 1 1 i 

1,060 

2 

0 

XlVa 

966 

• • 1 

2 

15 

XlVb 

680 

1 

0 

XlVg 

792 

2 

37 

XlVh 

754 

3 

151 

XIVi 

680 

1 

0 

XVg 

755 

1 

o 


NOTE: The mean field strengths have been rounded off to three significant 

figures. 
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The uncertainty in field strength corresponding to this wavelength 
uncertainty is ±75.5 gauss for }.5250.2. That value also represents 
the minimum field strength that can in principle be detected from these 
spectroscopic data. In practice, the minimum may be somewhat larger 
because the resolving power of the spectrograph and the distorting effect 
of the finite slit width may determine the minimum resolvable separation. 

The total probable error for the magnetic field values can be esti- 
mated as the root mean square of 180 gauss and 75. 5 gauss, or 195 gauss. 
The probable error in the spectroscopically determined field strengths, 
therefore, is approximately ±200 gauss. 

SUMMARY OF ERROR STUDY 

Magnetographs are, in general, subject to large systematic 
errors. The magnitudes and nature of the errors depend upon the 
particular type of magnetograph. Numerical error estimates rarely 
appear in the literature; errors in the range of 10 to 20 percent are 
generally considered to be moderate {Ref. 40). Most systems are con- 
sidered to be more reliable for field strengths less than about 1500 
gauss than for very large fields. 

It may be possible for a large systematic error to go undetected, 
although fortunately that is less likely now that results from many dif- 
ferent instruments can be compared. [In 1962, Stepanov, Shaposhnikova, 
and Petrova determined that Mt. Wilson measurements underestimated 
a field of 1000 gauss by no less than 700 gauss, on the average (Ref. 12, 
p. 205).] 

If observations in \5250.216 are used, all previous experience 
shows that the calibration must be done by comparing the signal levels 
to those measured for the same solar features in another line for which 
a calibration can be made. 
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Refined error estimates, of course, must wait until the system 
has been in use long enough for observers to collect a large mass of 
da^ta and to conduct careful tests for instrumental error. In making 
a preliminary estimate, no justification exists for assuming the field 
strengths will have accuracy better than 20 percent of the measured 
value. Therefore, ±20 percent is recommended as a provisional 
minimum error estimate. 
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VARIATION OF CIRCULAR POLARIZATION CONSTANT 
WITH POSITION ON THE SOLAR DISK 


For real-time rapid reduction of magnetograph signals, the 
relation between degree of circular polarization and magnetic field 
strength is assumed to be linear, for fields up to about 1000 gauss. 
The magnetograph calibration includes determining the slope for 
this relation. The value of Ci , however, depends explicitly upon 
position on the solar disk. 


If the Unno formulation for polarized residual intensities is 
used, the variation of Cj with position on the solar disk can be written 
as a function of known quantities (the Unno formulation is used because 
it permits C x to be written as a linear function of (3 0 , the limb darken- 
ing constant). The expected variation in Cj across the disk can then 
be estimated using a minimum number of line -dependent parameters. 


For small magnetic fields, the relation between the mean 
circular polarization, P v , and the magnetic splitting is assumed to be 
given by 


= (Vh cos i|<) 


(3) 


where 


Vj + 6 


dH(a, v) 


-W*o *1 0 


/ 


v i 6 


dv 


1 +’l 0 H(a. v)| 


T(v, vj) dv 


vj + 6 


/ 


1 + 


v; 8 


1 + t[ 0 H(a, v) _ 


T(v, v^) dv 


(4) 
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Equations (3) and (4) are taken from Reference 29, p. 1-17; 
Equation (4) is Equation (1. 14) of Reference 29. All symbols have the 
same meanings as in Reference 29 and will not be redefined here. 

Because v is independent of |jl and (3 0 , the three integrals in 
Equation (4) can be written in a form that is independent of position on 
the solar disk. 



and 


v. + 6 



T(v, v { ) dv 
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Then 


-Ko Vo A 
Ci s 

B + P 0 O D 


9 c t -j3 0 t| q A p 0 O 2 n 0 A D 

dp (B+p/3 0 D) (B + p £ 0 D) 2 


a Cj c, 

a n h 


0 O a c i 
B + p P Q D 


(5) 


B and D are easily evaluated, and p 0 is known. So Equation (5) can be 
used to find the change in Cj for different positions on the solar disk. 


Since the value of Cj at p = 1. 00 has not yet been established, 
the quantity to be calculated will be 


C, {p) 


Ci (p = 1 . 00 ) 


Approximately, 


C, <p) 

C, (p = 1 ) 


Q (p + Ap) /ac 


Ci (p = 1 ) 


_ 

9 JU 


*p 

(p = 1 ) 


( 6 ) 
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So if Cx (|jl = 1) is assumed to be known, 

Ci (k = 1)=C 1 (|i = 1) 

C, (|x = 0.9) _ C, ([L = 1) f 1 p 0 D 

Cl (p * 1) Cl (jx - 1) Cl (|X = 1 ) [o7? B + 0. 9 P 0 D. ^’° * °‘ 9) 

Ci Ih- * 0.8) Ci (n = 0. 9) Ci{fi = 0.9)[l P 0 D 

Ci (m- ~ i) “ c x (p = l) “ Ci (m. = i) [o7s ~ jT+'o. 8 p o D 

etc. 

To evaluate D, several spectral line parameters must be known: 

ax d - the Doppler width, which must be known in order 
to evaluate H(a, v) 

a - the damping parameter, which is also required 
to evaluate H(a, v) 

- the ratio of line -to -continuous absorption 
coefficients. 

None of these quantities is well known for X5250 in sunspots. 
Fortunately, however, the values found from Equation (5) are far more 
sensitive to the better-known quantities ^ and than to the three line 
profile parameters. The use of photospheric values of AX £>, a, and 
rj Q should not introduce serious errors. 

To investigate this point, a set of maximum and minimum 
values for Ci (p) was calculated by using the smallest and largest 
values of AXj}, a, and q which would be expected on the Sun, for any 


(0.9 - 0.8) 
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spectral line formed in. LTE, The results for X5250 should lie well 
within this range. The line profile parameters used for this calcula- 
tion are given in Table 5. The value of f3 0 was 1. 70 for both cases. 


TABLE 5. LINE PROFILE PARAMETER VALUES FOR 
DETERMINING MAXIMUM AND MINIMUM 
Ci (y)/Ci 0) 


PARAMETER 

VALUE FOR MAXIMUM 

VALUE FOR MINIMUM 

AXq 

■ " 

5.0 mini angstroms 

40.0 mil liang stroms 

i a 

0.05 

0.295 

% 

5.0 

30.0 


From Equations (5) and (6), it appears that the smallest possible values 
that could be obtained for Cj (p)/^ (1) would be found if p o = 0. A set 
of Ci (p)/Cj (1) values was calculated for this case; it is tabulated in 
Table 6 under the heading, "lower envelope". The maximum and min- 
imum values found using the parameters given in Table 5 also are given 
in Table 6. The difference between these sets of values is an over- 
estimate of the error that would be introduced by using photospheric 
parameters to calculate Cj (p) for sunspots. 

When p is less than 0. 5, the maximum error which would be 
introduced by using photospheric values is certainly only a few percent. 
It should be emphasized that the range of values covered in Table 5 is 
truly enormous, and the parameters for any single line can be estimated 
with reasonable certainty to within much smaller limits. 

The results in Table 6 suggest that the C* (p) values for sun- 
spot umbrae can be approximated by using the A X D> a » and rj 0 values 
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TABLE 6. RANGE OF VALUES OF C l (p)/^ (1) FOR PHYSICALLY 
POSSIBLE VALUES OF LINE PROFILE PARAMETERS 
AND FOR 3 0 = 0 


v 

MAXIMUM 

Ci (u)/c, (1) 

MINIMUM 
Ci (w)/C, (1) 

LOWER 

ENVELOPE 

j 1.000 

1.000 

1.000 

1.000 

! 0.917 

0.962 

0.929 

0.909 

0.833 

0.920 

0.855 

0.818 

0.750 

0.873 

0.778 

0.727 

0.667 

0.819 

0.698 

0.636 

0.584 

0.758 

0.614 

0.546 

0.500 

0.686 

0.526 

0,455 

0.417 

0.603 

0.433 

0.364 

0.333 

0.502 

0.335 

0.273 

0.250 

0.382 

0.232 

0.182 

0.167 

0.231 

0.122 

0.092 

0.083 

0.024 

0.003 

0.000 


which are found for X5250 in the photosphere. That approximation 
must be made because of the scarcity of high-quality umbral profiles 
of X5250 in unpolarized light. The value of (3 Q certainly changes 
between the photo sphere { and sunspot umbrae, and this factor is 
dominant. Fortunately, (3 0 for sunspot umbrae can be determined. 
Two separate calculations of Ci (h-)/Ci (1) were therefore made: 
one for the photosphere and penumbrae, and one for umbrae. 

For the photospheric calculation, the values of AXj-y a, and 
Tj 0 were found by fitting an absorption profile to M. J. Hagyard's 
observations of X5250 in the photosphere. The values determined 
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were AX^ = 12 milliang stroms, a = 0. 2, and tj 0 = 9. For the photo- 
sphere at X = 5250 angstroms, (3 q = 1.70 (Ref. 41, p. 171). According 
to Moe and Maltby (Ref. 42), M The intensity ratio between the penumbra 
and the photosphere shows little or no variation with position on the solar 
disk. " The same value of limb darkening constant (3 0 will therefore be 
good for both photosphere and penumbra, and the line profile parameters 
will be the same because good penumbral profiles of X5250 are as rare 
as good umbral profiles. The resulting values of Cj (}J-)/C| (1) are 
given in Table 7. 

TABLE 7. Ci (y)/Ci (1) FOR PHOTOSPHERE AND PENUMBRAE 


: y 

Cl (p)/C, (1) 

1.000 

1.000 

0.917 

0.957 

0.833 

0.908 

0.750 

0.855 

0.667 

0.796 

0.584 

0.729 

0.500 

0.652 | 

0.417 

0.566 j 

0.333 

0.464 

0.250 

0.345 

0.167 

0.202 

0.083 

0.016 


! For the umbral calculation, the same values of AXp, a, and 
T) will be used, but a new must be found. Wittmann and Schroter 
(Ref. 43) give center -to-limb intensities of sunspots at X 4680 
angstroms and 6400 angstroms. On the authority of Figure 3 of 
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Reference 44, the values for X = 5250 angstroms can be found by linear 
interpolation between the two sets of Wittmann-Schroter values. The 
results are shown in Table 8. 


TABLE 8. CENTER-TO-LIMB INTENSITY DISTRIBUTION FOR 
SUNSPOT UMBRAE AT X = 5250 ANGSTROMS 



I u (4680 angstroms) 

INTERPOLATED 
I u (5250 angstroms) 

I u (6400 angstroms) 

1.00 

3.68 

3.77 

3.94 

0.90 

3.66 

3.75 

3.93 

0.80 

3.63 

3.72 

3.90 

0.70 

3.60 

3.69 

3.87 

0.60 

3.56 

3.65 

3.83 

0.50 

3.52 

3.61 

3.78 

0.40 

3.47 

3.55 

3.70 

; 0.30 

3.35 

3.33 

3.59 

0.20 

2.98 

3.12 

3.40 

0.10 

2.24 

2.50 

3.02 


The change in intensity between center and limb is quite small. 
To obtain (3 Q , the values of I u (p)/I u (p = 1) are fitted to the 

form 

i u M 

=1 • U] + (Uj) (n) 

I u (D 

where 

u i . . . ■ ■ 
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A simple equation of the form of Equation (7) does not permit a good 
fit over the entire range of p values, but an acceptable fit for 1. 0 ^ p 
^ 0.4 can be obtained by setting p o = 0.093 (uj = 0.086). In Table 9, the 
values of I u {p)/I u (1) for p o = 0. 093 are compared to the interpolated 
measured values. 


TABLE 9. VALUES OF UMBRAL LIMB DARKENING FOR 
3 0 = 0.093 COMPARED TO MEASURED VALUES 


1 

1 y 

iu (y)/iu 0 ) » 
3 0 *; 0.093 

Iu (|0/I u 0) 
FROM TABLE 8 

1.00 

1.000 

1.000 

1 0.90 

i 

0.991 

0.993 

0.80 

0.983 

0.987 

0.70 

0.974 

0.978 

0.60 

0.966 

....... ■■ 

0.968 

0.50 

0.957 

0.958 

0.40 

0.948 

0.942 

1 0.30 

! 

0.940 

0.884 1 

i 

• ! 

[ 0.20 

0.931 

0.828 

1 

1 0.10 

0.923 

0.664 

I 

< 

[3 value is less accurate 
0 

in this range) j 

1 


The value (3 0 = 0. 093 was used to calculate Ci (p)/Cx (1) for 
umbrae. The values of Cj when p is less than 0.40 should be given 
less weight because the approximate limb darkening relation is less 
accurate in that range. 

The line profile parameter whose value in umbrae is most 
uncertain is probably r| 0 . Values Cj were calculated using t) q = 9, 
the photospheric value, and ri 0 = 40, since tj q may be much larger in 
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umbrae. The results for both r jo values are given in Table 10; they 
are almost identical. 


TABLE 10. Ci (y)/Cj (1) FOR SUNSPOT UMBRAE 


y 

Cl (u)/Ci (1), no - 9 

Ci (y)/Ci (1), n 0 = 40 

1.000 

1.000 

1.000 

0.917 

0.915 

0.913 

0.833 

0.827 

0.825 

0.750 

0.740 

0.737 

0.667 

0.652 

0.648 

0.584 

0.562 

0.558 

0.500 

0.471 

0.467 

0.417 

0.380 

0.376 

0.333 

0.286 

0.283 i 

■ 

0.250 

0.192 

0.190 


0.167 

0.098 

0.096 


0.083 

0.000 

0.000 



(3 Q value is less 

accurate in this range) , 

8 


Should it become necessary to do so, a supplement to Table 10 
for fa. < 0.4 can be calculated, using a value of (3 Q which gives a better 
fit to the umbral intensities near the limb. 

The values of Tables 7 and 10 are shown plotted in Figure 5. 
The relation for the umbra is almost linear. 
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PENUMBRA 





PROGRAM TO LOCATE AND PLOT MAGNETIC 
FIELD NEUTRAL LINES 

Two of the properties which are of the greatest interest in solar 
magnetic field maps are the regions where the magnetic field gradients 
are steepest and the locations of the boundaries between field regions 
of opposite polarity. The regions of steep gradient are easily located 
on magnetic field contour plots because they are the regions where the 
contours are most closely spaced. However, the boundaries between 
opposite -polarity field regions (the boundaries will be called "neutral 
lines" for brevity) are generally less distinctive in appearance on 
contour plots, and it is convenient to provide a separate program to 
plot the neutral line locations directly from magnetograph data. 

A FORTRAN program to make such plots has been written for 
the Univac 1108. The position of each neutral line is identified by 
selecting the coordinates of pairs of adjacent points for which the 
measured signals are opposite in sign. The program in its present 
form will analyze a square grid of data 128 points on a side and will 
store 1998 coordinate values. Places where the signal changes sign 
because of the sign reversal of the magnetograph calibration curve for 
large magnetic field values will be identified by the program as portions 
of neutral lines. When the field strength values at which the calibration 
curves change sign have been established, the program will be modified 
to ignore signal sign changes for field strengths in the reversal range. 

The program compares points in two adjacent horizontal data 
rows. Beginning with rows 1 and 2, the signs of pairs of points are 
compared by taking groups of four data points at a time. The points 
are labeled as shown in Figure 6; their signs are compared in the 
order PI, P2; PI, P3; P2, P3; PI, P4. If all four points are of the 
same sign, 
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nw9i e. mm uselihg patterns for 
nm?.PL like mmm 

‘{/Ant# H &fld F4 &re relabeled FI and PE, and the next vertical pair 
to the right 3P3 and P4 for a new four-point group. The 

',or/)pti,rt&Oft proc&nn in repeated across the data row. When two points 
of opposite £%n are found, their coordinates are stored. In order to 
#aye n phots in. memory, the program then stops testing values in a 
particular four -point group and moves to the next group. This explains 
v/hy, in the final plots, thf? lines which connect points of opposite sign 
may i&he any of the four directions; they simply record the relative posi- 
tion of the members of the first pair of points of opposite sign which was 
found in each group Of four points. When the first two data rows have been 
scanned, the process is repeated with the second and third rows, and the 
full MQqVkQnGQ is repeated until the entire frame has been scanned or until 
the array reserved for coordinate values has been filled. 

After the scanning has been completed, the coordinate values 
of the pairs of points of opposite sign are plotted by connecting them 
With short 'lines. A lino drawn joining the midpoints of these short 
eomtecting lines will trace the location of a neutral line to the precision 
allowed by the spatial resolution of the data. To assist in determining 
the coordinates of the plotted points, grid lines are drawn at every 
fourth data row and column. 
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The program was tested on a frame of magnetograph data which 
was obtained from the Kitt Peak National Observatory and which has 
been used for testing other magnetograph analysis programs at MSFC, 
This particular data set forms a grid 120-by-120 data points in size. 
Three cards were missing from the data deck; one of these contained 
data which crossed a neutral line and produced a gap in the upper 
right-hand corner of the final plot. Except for this gap, the plot agrees 
in every detail with the neutral line pattern obtained from the data by 
the Kitt Peak Observatory. It also agrees with a hand plot of the data. 

The plot of the test data is given in Figure 7. Figure 8 is a 
flowchart of the test program, which reads data from cards, and Figure 
9 is the listing of the test program. For use with real MSFC magneto - 
graph data, the program will be modified to read the appropriate 
magnetic tape formats. 
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FIGURE 7. PLOT OF NEUTRAL LINES PRODUCED BY TEST PROGRAM 
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FIGURE 8. FLOWCHART OF TEST PROGRAM 
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C o?£PC GPICIFNT 5 CONTOURS IN HI £ NF T OP R IPH till 

C VEP'iCN FCP TEST WITH KPNO MIGNETOGPIPH CITI 

CIMENSION > I ( 1 2 H > , X ? ( 1 ? 9 ) • NLXI?OOfll, NLYI 2 UOC) 

Cl MENSI CN JPRIV 122 ) 

C»TI IIRRI VI I ) #1 = 1 * 221 / 2 MHHIRCC 0 PV CNL Y * ONE COPT « 18 *EP 
t<PE= 5 
KWP= G 

RE 1C I K P E • 11 RMSN 

C REIC RMS NOISE LEVEL FOR FR * ME 

1 FORMlT I F IC.C) 

ILCCr 1 

C CCUNTEP FCR SZERCff GR *E I EN T S 

REIC IK p £ * 051 IMIX. IM*Y. ISTZF 
M 5 FCRM 4 T 13110 

C NC. OF ROWS* NO. OF COLUMNS* NO. OF NF U T P* L LINE POSITIONS TO 8 F 

C S TC P EC 

IS 21 = I p I ?E - 1 

C REIC 1ST 2 p CWS CF CITI 

REIC IK C E * MG) 1X1 If). I - 1. I Ml X I 
PE I C IKPE* ME) I X 2 I I ) • ir 1, IMI X I 
U£ FCPMIT | J2IF5.0* 1X1 1 

C SET Y CCCPCIMTE COUNTER INC STIR! SEIRCH OF cpl^F 

NM|V= I M I Y * 1 

CP 9 NYI: 1 • NMIY 
N Y 2 - NYI *1 

C STIPT CF iGPICIENTa EVILUMION LOOP 

NMIXr I M l x - 1 
CC IP I = 1* KMIX 

C SET X CGCRCTNITE COUNTER 

NX 1 = I 
NX 2 = 1 ♦ 1 
PI r XI f I ) 
p 2 : X 2 1 1 1 
Pi: XIII *1 1 
PR: X 2 1 1 ♦ 11 
IF IPU 1 CM • 3 . 1 C? 

101 IF 1 P 2 1 M* ?. 3 
10 M IF I P? 1 3 * ’, M 
M NLXIILCClz NX] 

NLYIILOC): NYI 
i'LXULOC ♦ Tl: NXI 
NL Y II L 0 C •* 1 1 = NY 2 

WRITE IKWR, 301 PI, P 2 » N XT, NYI 
7 L F OPM IT f lH , 2 1 E 1 2 • * » 5 X 1 , 2 IIP, 5 X 11 
I L CC: ILOC *2 
IF ULCC-IS 21 ) IP, e , 5 
5 IF INM -701 E, 7 » 7 
F WRZ TF I KWR, B) NYI 

B FPRMIT I 1 H • 59 HIRRIY R E C FH V£ C FOR 7 ER 0 GRIC 1 FNT COORCINUFS IS Ft 
ltLEC NYI:, 15 ) 

WRITE IKUR* 111 

It FCR*IT I 1 H , IIHTERHINITE RUN) 

GC TO 9«9 

I IF IP I 1 1 12* 12 * I) ? 

112 IF IP’) 12 * 12 * II 

FIGURE 9. LISTING OF TEST PROGRAM 
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11 ? 
I ? 


1? 
115 
114 
1 4 


15 

IIP 

I1C 

If 
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I i 6 
1 1 7 
1 7 


IF IP3) 1 ? • 12. 1 2 
NLXtlLOC): N X 1 
NLY II LCC): KYI 
NLXI ILOC ♦ 1 ): NX? 

NLY II ICC ♦ 11= NY 1 

WRITE I K UR * 30) PI. P T . NX1# NY1 

ILCC: ILCC « 2 

IF UL 0C-IS71 ) 10. e . 5 

IF IP 2 ) 115. 15. 114 

IF IP*) IS. 15. 1« 


115 

1.06 


18 
1 5 

?n 


7 

172 


IF IP?) 14* 1 


IF 


NL XI It OC ): N X 1 
NLY II LPT ) = NY 2 
NL XI IL 0 C ♦ 1): NX? 

NL Y II LCC ♦ 1 1= NY1 

WRITE IKWR. 30) P2, P7, NX1. NY 1 

I LCC : ILCC ♦ 2 

IF I IL CC-IS71 1 IP. c . 5 

IF IP1) l lp 9 IP. IIP 

IF < PH 1 16 * 10. IP 

IF f P 4 ) 1C. IP. IF 

NLXI ILCC )- N X 1 

NLY II LCC I = NY1 

NLXIILOC ♦ 1)= NX? 

NLY IILCC ♦ 1 )= NY 2 

WRITE CKWR. 70) PI. P4 . NX1. NY 1 

ILCC = ILCC ♦ 2 

IF 1 1LCC-IS71 I IP. 5. £ 

CPNTINUE 

PI: XI IIMIX) | 

F2 : X ? ( I M 4 X ) 

IF I PI) 119. 18. 117 
IF IP? ) 18. 18 . 17 
IF IP?) 17. 18 . 18 
N L X IILCC): I FIX 
NL VI ILCC )= NY1 
NL X II LCC ♦ 11= IMIX 
NL Yl ILCC 4 1 )= N Y2 
ILCC : ILCC 4 ? 

IF I 1LCC-IS?! ) 18. 1 8. IP* 

IF INY 1-70 J 1PF. 7. 7 
WRITE IKWR, 8) NY1 
WRITE (KUP. 11) 
cc TC 999 

IF INY1-NK4Y) 1 * . 7. 7 
CC 20 12= 1. 14*4 X 
XI III): X?l I I 7 
C CNTINUE 

RE 4 C NEXT RCW CF Cl T4 

RE4C |KR£, 4 f ) (X2CI). I: 1. IMIX) 

CCNTINUE 

T HI * CCHPLETF5 5EIRCH FOR S 2€ RO C R4 CIFN? S POINTS 
IF IN: ILCC - I 
WRITE IKWR. 132) I FI N 
FORMIT 1 1 H . ShIFIN:. 

WRITE IK WR . 21 1 


IIP ) 


original eage^ 

OP POOR QUALHi 


FIGURE 9 - Continued 
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21 FORWIT UH .35KX CCCRC OF 2ERO CRICTENT POINTS I RE 1 
WRITE <K WR t 27 ) CNLXlDt 1= It IF I N) 

77 FCPPIT ( 2 H t 54 I5t 5X ) » 

WRITE (KWRt 23) 

23 FOPHIT CIH t 3SHY CCCRC OF ZERO GRICIENT POINTS I Re ) 

WRITE 4KWR# 2?) (NLY(I)t I? It IFIN) 

PLOT TIPE 

C 4LL KENT 4 9 1 IRRIY) 

cut 6 Ric.iv tr tu.ii.i ?o.Ptf .u *i2r.r. m>*p.r *-5*5*- ic»-ic #♦*•♦«» » 

I Fit ILCC-2 

CO BOO 1= It I F 1 1 2 

X: j NL X f I ) 


ciU 

XSCLVUXt 

NX P t 

I ERR I 

NX I t 

NXP 



x t , Nt x 1 1 + n 



C*tl 

xfclvh Xt 

NXP t 

1ERR) 

NX 2 1 

NXP 



Yt 4NN«Y«U-NL Y(D 


CALL 

YSCLVK Yt 

NY P t 

1ERR > 

NYlt 

NYP 



Yt INHIY* 1I-NIV4 1* 1 

1 

cm 

Y*CLV 1 4 Y t 

NYP t 

IERR) 

NY 2 1 

NYP 




CHI L INEVlNXltNYItN X?tNY2) 
gnO CONTINUE 

CIU ENC4CB 
929 5 TC e 
ENC 


FIGURE 9 - Concluded 
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CONCLUSIONS AND RECOMMENDATIONS 

Even the best magnetograph results are subject to large uncer- 
tainties. It is unwise to assume the uncertainties in results obtained 
with a new system are less than ±20 percent. They may be larger. 

To reduce some of the uncertainties of calibration and inter- 
pretation, it is recommended that the magnetograph be modified to 
operate using a line other than the neutral iron line at 5250. 216 
angstroms. The new line, of course, should be carefully selected to 
introduce a minimum number of new uncertainties. If compelling 
instrumental reasons exist for requiring the new line to be near 
5000 angstroms in wavelength, both A5253 and X5324 would be worthy 
of closer investigation. Because the amount of Zeeman splitting is 
proportional to the square of the wavelength, lines of wavelength 
longer than 5000 angstroms would tend to produce larger splittings and 
therefore could be advantageous choices. A number of red spectral 
lines have been suggested by various authors as possible magnetograph 
lines; two neutral iron lines which should be investigated are A.6173.348 
(lower excitation potential = 2.213 electron volts) and X6302. 508 (lower 
excitation potential * 3,671 electron volts) (Ref. 12, p. 179). A number of 
other candidates may exist. 

Evaluation of the Characteristics of possible new magnetograph 
lines should begin with a preliminary study of the type summarized in 
Table 3. Lines which are not rejected in the preliminary study must 
then be further examined to prove that their properties arc compatible 
with the assumptions implicit in the magnetograph calibration analysis. 

As was shown earlier in this report, attempts to calculate line profiles 
for sunspot models may reveal such incompatibility. The profile 
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calculation thus forms a test which should be applied to each line which 
is seriously considered for magnetograph operation. 

When the magnetograph is operated using X5250, it is recommended 
that no reliance be placed on calibration systems which require, even 
indirectly, the calculation of residual intensities of that spectral line 
in magnetic fields. 
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SECTION 3 

PROGRAMS FOR PRESENTATION OF 
VECTOR MAGNETOGRAPH DATA 
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INTRODUCTION 


This section contains documentation for four programs written 
for the Univac 1108. Two of them have been described in Teledyne 
Brown Engineering Interim Report No. EE-MSFC-1815 (May 1974). 
Some minor changes have been made to these programs, so updated 
descriptions and current listings are given here. These two programs 
plot selected portions of Real Time Solar Magnetograph data: one 
program displays locations of magnetic polarity changes, and the other 
sHows relative steepnesses of magnetic field gradients as a function of 
location in the data matrix. The third program removes isolated large 
data values from the data array and replaces each of the removed 
values by the average of the values surrounding it. The "corrected" 
array is then written on a new tape in the same format as the original 
data. The fourth program is a simple routine which combines the "U" 
and "R" matrices from the original magi^tograph data to form the 
matrix of values corresponding to transverse magnetic field data. 

This new matrix is written on a new tape in the same format used for 
the longitudinal data. 
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MAGNETIC FIELD NEUTRAL LINE PROGRAM 

This program locates and pu cs the locations in the magnetogram 
at which the magnetic field changes polarity. At these points, the sign 
of the magnetograph signal changes, so the program really plots the 
locations where the signs of adjacent signal values are different. 

The position of each neutral line is identified by selecting the 
coordinates of pairs of adjacent points for which the measured signals 
are opposite in sign. The program in its present form will analyze a 
square grid of data 128 points on a side and will store 1998 coordinate ; 
values. Places where the signal changes sign because of the sign 
reversal of the magnetograph calibration curve for large magnetic 
field values will be identified by the program as portions of neutral 
lines. When the field strength values at which the calibration curves 
change sign have been established, the program will be modified to 
ignore signal sign changes for field strengths in the reversal range. 

The program compares points in two adjacent horizontal data 
rows. Beginning with rows 1 and 2, the signs of pairs of points are 
compared by taking groups of four data points at a time. The points 
are labeled as shown in Figure 1; their signs are compared in the order 
PI, P2; PI, P3; P2, P3; PI, P4. If all four points are of the same 
sign, points P3 and P4 are relabeled PI and P2, and the next vertical 
pair to the right becomes P3 and P4 for a new four-point group. The 
comparison process is repeated across the data row. When two points 
of opposite sign are found, their coordinates are stored. In order to 
save space in memory, the program then stops testing values in a 
particular four -point group and moves to the next group. This explains 
why, in the final plots, the lines which connect points of opposite sign 
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data row 1 

• 

• 

PI 

. . . • 

P3 

• 

• 

data row 2 





• 




P2 

P4 



FIGURE 1. POINT LABELING PATTERNS FOR 
NEUTRAL LINE PROGRAM 


may take any of the four directions: they simply record the relative 
positions of the members of the first pair of points of opposite sign 
which was found in each group of four points. When the first two data 
rows have been scanned, the process is repeated with the second and 
third rows, and the full sequence is repeated until the entire frame has 
been scanned or until the array reserved for coordinate values has been 
filled. 

After the scanning has been completed, the coordinate values 
of the pairs of points of opposite sign are plotted by connecting them 
with short lines. A line drawn joining the midpoints of these short 
connecting lines will trace the location of a neutral line to the precision 
allowed by the spatial resolution of the data. To assist in determining 
the coordinates of the plotted points, grid lines are drawn at every 
fourth data row and column. 

The program uses as data the magnetograph output tapes in 
their format after undergoing initial processing. In this format, the 
transverse and longitudinal components are written on separate tapes. 
Each magnetogram consists of 9 records: a label record of 19 words, 
and then 8 records of 2049 words each (a counter plus 2048 signal 
values). In some cases, a few words at the beginning of each mag- 
netogram were lost when the original tapes were processed. The tape 
read portion of this program contains a short routine to correct for 
this condition. 
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r Two values must be supplied in order to execute the program. 
Both of these are entered in DATA statements within the program. 

Th£ first is the number of longitudinal matrices contained on the tape. 
This quantity is given the variable name of IPP and is entered as an 
integer variable. For example, for a tape containing 12 longitudinal 
matrices, the DATA statement would have the form 

DATA IPP/12 / 

Only longitudinal matrices are processed by this program. The determi 
nation of transverse field strengths from data recorded in the format 
used here requires that a transverse field value be calculated as the 
square root of the sum of the squares of its components. The algebraic 
signs of the transverse field values therefore cannot be determined, 
and polarity reversals cannot be detected. 

The second datum to be supplied is the mean background noise 
level for the frame of data. This variable is called RMSN; for the 
matrices that have been processed so far, the value RMSN = 0.05 is 
usually satisfactory. RMSN is a real variable. The DATA statement 
which specifies RMSN would be written 

DATA RMSN/. 05/ 
for a typical data frame. 

A listing of the program is given in Figure 2. A sample plot 
is shown in Figure 4. 
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WTira 

1* 

c 


• 2EB0 SR AOrcWT* - 'COTTTOOWT In 3TTSGWto6TFAPH DAT* Mofl J/74 ' ““ 

00 tot 

2* 



DIMENSION Ni.7<2000> tNLVl?0nol •ARRAV122) f TEMP(6000l tY<32, 1281 

00103 

3* 



DIMENSION Lt1ME(9) 

00101 

1* 



0IMEN510N A ( 20 ) 

OC105 

*• 



0)|ta ( array < I 1 , I *1 , 22 )/2lHH*R0C0PY ONLY, ONE 00*7 ,IB*6H / 

00107 

6* 



DATA I PP / | 2 / 

001 1 1 

7 • 



D/iTA HU / * / 

00113 

8* 



DATA RMSN/.nS/ 

001 is 

9* 



KRE- S 

OCI 1* 

10* 



«*R* 6 

00117 

II* 



CALU OPEN ( p , l , S> 

-ggi?0- 

12* 



IMAX. 128 

00121 

13* 



1 M A Y* 12« 

00122 

11* 



ISI2F- 200 ? 

00122 

15* 

c 


NO* OE ROWS, NO, OF COLUMNS, NO, of neutral LIN E POSITIONS to be 

00122 

1 6* 

c 


STORED 

00123 

17* 



I SZ I • 1SI7E - 1 

0012** 

1 ft* 



CAUL IDFNj ( 9, ARRAY) 

00121 

!<»• 

c 


OATa ape WRITTEN ON TAPE FROM right to LEFT, IN RON OROFR, ft 

0CI7H 

20* 

c 


rows at a time, precfoeo by a record containing 19 noros 

001 2H 

21* 

c 


OF HFAOER INFORMATION 

0C12S 

22* 



00 999 IP. 1, Ipp 

001 25 

23* 

c 


READ HEADER INFORMATION 

00130 

2 H • 


3*0 

CALL RF.DTPR 18, 2, |ER, N*, 19 , A) 

00131 

2** 



IF CMW-19) joo, t, 300 

00131 

26* 


1 

WRITE IKWR, 9901 NW 

00137 

27* 


990 

FORMAT 1 1 H , 3 H N w « , 151 

00137 

28* 

c 


RMSN is H MS NOISE LEVEL FOR FRAMF 

00110 

29* 



IF 11ER - 1 > 3 9 1 , 390 , 3’| 

00113 

30* 


39 1 

WRITE IKWR, 3921 | FR , NR 

00117 

31* 


3*2 

FORMAT (1H , 2 1 1 0 1 

00117 

32* 

c 


read |ST 2 rfcoros pr maGnftogram 

OP ISO 

33* 


39fl 

DO 65 !■ 1 , 2 

001*3 

39* 



ICT . 1 

001*1 

35* 



write 1KWR, 9881 ICT 

001*7 

36* 


9®ft 

FORMAT MH , 9H1CT-, IS) 

00160 

37* 



J. 2018* ( I-) 1 *1 , ~ 

OCI 61 

38* 


1 9 

CALL RFOTpr <8,2, IER, NW, |, 1 CNT , 2090 , TFMP 1 J > 1 

00162 

39* 



WRITE ( K W R , 98 9 ) NW, |fNT 

OCI 66 

10* 


9 69 

FORMAT (1H ,3HNW«,17,S>-,5HfCNT»,t*) 

00167 

11 * 


s 

IlF tNW-?"«19) I3r, 129, | Jo 

00172 

1 ?* 


1 30 

WRtTF ( KWr , 3921 NW, IfNT 

00176 

13* 


| 

K'K* 2019-NW ♦ 1 

00177 

11 * 



JX* 1 

00200 

IS* 


} • 

DO 1 2 1 *■ Kk, 2019 

0020 3 

16* 



t;emp ( K » * TFMP( JX 1 

00201 

*7* 



JX » JX ♦ 1 

00205 

1 ft* 


121 

continue 

Q0207 

19* 



* 2 - KK * 1 

002)0 

SO* 



00 1 3| X* 1 , K 2 

002 I 3 

Si* 



TEMPIKI • 0 . 
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00218 

52* 


131 

CONTINUE 



002 1 6 

53* 


li29 

IF CIER-11 395, 8650. 195 

......... 


00221 

SH* 


395 

WRITE (KWR, 392 1 IER, NR 



0022 $ 

55* 


8 65q 

IF <1-1 > 65, 865, 65 

. 


00230 

56* 


865 

IX • 100 



00231 

57* 




I Y* IT - 5c ” 


— 

00232 

56* 




IDELY- 25 



00233 

5?* 




N« FLO ( 2 , Ip, TFMP1 1 1 1 



00238 

60 * 




! 0 AYS* 200*FLD(26,| ,N 1 ♦ | 00*FLD I 2 7 , j f N J + a O*FLO 1 2« » 1 ,N> 


002 3 8 

6 1 • 




I ♦83»FL0(29, f ,N»* 2 o*FLD 130 , | ;m>*10*FL«!3 1 , 1 , N 1 ♦ 8 »Fl 0 1 32 . l . M » 


0023*4 

62 * 




2*8*FLD<33 , 1 ,NMp*FL0(38, I ,Nj*FLD«35, 1 , N ) 



00235 

63* ' 




N« FLO • 18,2, TEMPI 1 11 * - 



002 3 6 

68 * 




FLD<3p,8,Nl» FtDl I2,8,TFHP( | 1) 



00237 

65* 




1 HOURS* 2p»FLD< 30, 1 , W » ♦ 1 0*Fi.D t 3 I , 1 ,N.Wn«FLO| 32 

i | ,N > 


00 23 7 

66 * 




l*8*FLD(3 3,| ,N|*?*pL0t3«, 1 ,M1*FLD«35 ,I ,n» 



002 HO 

67* 




N* FLD1 20, 7 , TEMPI 1 ) > 



00 2 H ] 

68 * 




I H IN* 80*FL0<29, 1 ,N1*20«FLO(30: l , N 1 ♦ 1 0*fL D ( 3 1 , 

1 ,N> 


0C2H1 

69* 




I ♦ 8 *FLI>I 32, 1 ,N>*8*FL0l33 , I , M'» ♦ 2 *FrDT3T, • N ) ♦ f 1 FTT5TT7F5 


CC 2 H 2 

70* 




N* FLO ( 27 ,7 , TEMPI 1 ) ) 



002H3 

71* 




! SEC* 80*FLp«29, I ,N»t2p*FL0( 30, 1 • N > ♦ 1 0»FL D ( 3 l , 

I , N 1 ♦ 8 • F L 0 1 1 2 , 1 ,N» 


002H3 

72* 




l* 8 *FLD (33, | ,M)*?*FL0« 38, 1 ,N)*FL0»35 , 1 ,N> 



0P2HH 

73* 




ENCOHt ( 80 , L T I ME 1 1DAYS, INCURS, IMlN, ISEC 



00252 

7 H • 



HQ 

FORMAT (18, 5H PAYS, 17, 6 H HOURS , I 8 , 8 H M , NUT E 8 , 

18, 8H SFCOMos* 


0025 3 

75* 




CALL FRA M EV|3) ~ 



00258 

76* 




CALL PRINTV (58, LT 1 ME , lx, !Y 1 



00255 

77* 




call stoptv 



00256 

76* 



6 ft 

continue 



0024.0 

79* 




1L0C* 1 



00261 

* 0 * 




IC* 1 



00761 

81* 

C 



!CT IS NO, OF RFCORDS REAP' FROK TaPT“ 

— 


00261 

82* 

C 



tc is mo. of times read from Tape is executed 



00261 

83* 

c 



IC SM 0 UL n 8 F, EQUAL TO fcT . 1 



00262 

89 * 




f • 129 



00263 

85* 


.. i 


00 68 jj* p, 3 ? 



00266 

86 * 




00 70 fK* 1 , |20 



0027 1 

87 • 


1 


Y ( J J , KK » • TEMPI J ) 

„ ■ 

— 

00272 

88 • 


| 


1 * ! ♦ 1 



00273 

89* 



70 

CONTINUE 



00275 

90* 


|68 

CONTINUE 



00275 

9 1 • 

c 

j 


remove STRtpF 



00277 

92* 


j 


no 8P56 I M* 1 , 32 



TO3PT" 

93* 




IX* (I M-l ) *| 28 ♦ 1 |1 — 



— 

00303 

9<*« 


! 


TEMPI | x 1 * TEMPI 1 X ♦ | | 



0P3pM 

95* 


8056 

continue 



00 3 0 H 

56* 

c 

! 


SET 1 X COORDiMaTf COUNTFR AMp START SEARCH OF FRAME 


00306 

9 7* 




N M A X * T M AX - 1 



00307 

98* 


: 


11 * 2 



“ '0*307 

99* 

c 



start searching columns p tmROmgh i 6 



— 

00310 

100 * 


1 


N M A Y • 1 M * Y -| 



00310 

lot* 

c 



SPIKE NOISE FILTER 



003) I 

102 * 


21 1 

no 8060 JJ* 1 , 31 



003 |H 

103* 




00 8051 KK* 1 , 1 27 



03317 

108* 




X TEM* A8S(PmSN-A85( y( JJ,KK 1 ) 1 



00320" 

1 05* 




IF f»TFM-3,,RHSN> 8 o5 1 , *»05|, t<o51 ~ — 

— 

• 

00323 

106 * 

805 2 

YTEM* AflSIRMSN-ARS! Y( JJ*l,Kx» ) 1 



CO 32 41 

1 07 • 




IF I YTFH-3..RM5M) 8pS3, SqSj, 8q5i 
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00327 

108* 

00330 

109* 

0033 1 

110* 

00332 

III* 

00335 

112* 

00336 

fl3* 

003 *4 1 

1 19* 

003*43 

115* 

00385 

I 16* 

0,03*45 

1 17* 

00385 

I IS* 

003*45 

1 19* 

00350 

120* 

00351 

l ? I • 

0035? 

1 ?2* 

00352 

123* 

00355 

L?*4* 

00356 

1 25* 

00357 

1 26* 

00360 

1 27* 

00361 

128* 

0036? 

1 29* 

— 00 3-6 3. 

.1 30* 

00369 

1 3 » • 

C|0367 

132* 

00370 

1 33* 

00373 

139* 

00379 

I 35* 

003.7,5 , 

136* 

00376 

I 37* 

00377 

138* 

00802 

1 39* 

00903 

190* 

00906 

t 9 | • 

00807 

1 92* 

009| 6 

1 9 3 • 

009 | ! 

| 9*4* 

00*4 1 2 

1 95* 

008 15 

1 96* 

009»6 

1 97* 

00971 

198* 

00822 

199* 

00823 

1 50 • 

00929 

154* 

009?5 

152* 

00930 

153* 

00931 

159* 

60*439 

155* 

00935 

156* 

009 3 6 

457* 

00937 

4 58* 

00992 

159* 

009*45 

1 60* 

009 50 

16 1 • 

00951 

162* 

00952 

163* 


*<05*4 Y(JJ,KK) ■ RMSN/!.?o 
GO TO 9051 

•4053 ZTEM* AB5(RmSN-a8s< Y 4 >♦ > 

If I ZTEM-3,*RMSN> *4055, 9os 5 , r 0 5| 

9P5S 7 Z T E M ■ ABS 4 rMsN.AbS (YJ JJ *1 »**♦! ) » > 

If 4ZZTEM-3.*RMSN j <405*4 1 9059, "*405”i ~ ~ ” 

<4051 CONTINUE 
•4050 CONTINUE 

00 9 J J* 1J,|6 

C PREPARE To gET correct coordinates values TO BE stored 

c M i sunoerst ano i no about order in *hi c u tape has written 

C SO MUST Exchange X AND Y COORO AFTER ImE* ARE OFTERhINED 

NX 1 * JJ + 4 JC - I ) *16 
NX2* NX! ♦ j 
00 10 I • |, 427 

C l SET Y COORDtNATF COUNTER 
NY I ■ I 

I NY2* 1 ♦ 1 ' 

! PI" Y4JJ, I) 

P2- Y 4 J J » | ♦ | ) 

P3» Y ( J J* J , I) 

P‘4* Y ( JJ* | , 1*| ) 

PP* AB54PJJ-RM5N 

If 4 PP ) P02, 90? * RD3 - " ~~ 

p 02 PP* C, 

8Q3 IE 4 p I ) 9o<4 , 90?. 903 
909 P|« -pp 

GO TO 90? 

903 PI* PR 

902 PP* ABS1P21-RMSN " ~ 

IE 4 PP J Bo<4 , 901 . 805 
80*4 PP* O. 

805 IF 4P2I 9 0 5, 901. 9|6 
905 P 2* -PP 

GO TO 901 

916 P2* PP 

901 PP* ABS 4 P3 » -RMSM 

IE 4 PP I 806, 912, 807 

806 PP* 0. 

807 IE 4 P 3 ) 9 | 7 , 91?, 9q8 

917 P3* -PP 

GO TO 9 , 2 -• 

908 P 3* PP 

912 PP» ABS < P<4 1.RMSN 

IF 4 PP ) 808, 91 1 . Ap9 
BOB PP* 0. ::: ’“ 

«09 IE (P‘4) 909 , 9 1 1 , 9 |C 

909 pr* -pp •••'“ — 

GO TO 911 

910 PR* Pp 

91 I IF 4 P 1 ) IC*4, 3, IC3 
IP3 if 4P2) ,, 3 

109 IF 4P2) 3, 3, *4 

9 NLX4TLOO* flYl ... — - 

NL'Y.I ILOC 1* NX 1 
NtX 4 IL°C ♦ | 1 • MY2 


FIGURE 2 - Continued 
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OO **53 

16*4* 


NLY( ILOC ♦ 1 )- NXi 

- 00 **S** 

165* 


WRITE (KWR, 30) Pi, P2, NXI, NY1 

00 ** 6 2 

166 * 

30 

FORMAT «1H , 2 ( F I 2 • 5 , 5 X ) , 2 <! 6 i 5X1) 

GO **6 3' 

167* 


ILOC- ILOC ♦ 2 

00** A** 

166* 


IF (ILOC - fSZl) 10, 5, 5 

" 00** A 7 

' YTTW 

5' 

IF rNVI-70) 6 , 6 t 966 

00** 7 2 

1 70* 

6 

WRITE (KWR, 8 ) wTl 

00** 7 5 

! 7 ! • 

8 

FORMAT ( 6 H NY!** I 7 * 

00** 7 6 

172* 

1 

WRITE (KWR, *411) 

00500 

17 3* 


WRITE (KWR, *412) 

00502 

1 7*4* 

*411 

FORMAT (2SH ARRAY RESERVED IS FILLED) 

~ ' 00 5 0 3 _ 

|7S* 

*412 

Format KI5h tErMinaTS rUn> 

0050 M 

176* 

906 

IF (ICT- 8 ) 907. 7, 7 

“ 005 0** 

177* 

C 

60 TO END Or MAGNETOGRAM 

0050 7 

178* 

907 

CALL REDTPR (8,2*IER,Nw,i ,!CMT , 20*49 , TEmr( i ) ) 

00510 

1 79* 


WRITE (KWR, 989) NW, ICNT 

005 1 *♦ 

180 * 


ICT - !CT ♦ 1 

60S ! 5 

J 81 * 


WRITE (KWr, 988 ) ICT ~ 

005 20 

182 * 


TF f I C T - 8 ) 907* 7, 7 

005 20 

183* 

C 

GO TO END Op MAGWeToGRAM. THEN START Nr X T ONF 

005 2 3 

1 8 ** • 

3 

IF (PM 1 12, 12, 1 13 

005 2 6 

1 95* 

1 12 

IF (P3) 12, 12, 13 

005 3 1 

1 96* 

113 

TF (P3> 13. 12. 12 

0053** 

1 A7* 

! 3 

nlx ( tloc ) « n'Y i 

00535 

198* 


M L,Y ( 1 LOC ) * MXl 

005 3 6 

189* 


NLX (ILOC ♦ i )■ mV i 

00537 

i’C* 


NLY( ILOC ♦ i )- M X 2 

005 <4 0 

191 • 


WRITE (KWR, 30) Pi, P 3 , NX), NY 1 

005*46 

192* 


ILOC* ILOC ♦ 2 

605*47 

1^3* 


IF ( tLOC-lSz I ) 10, 5, 5 

00552 

1 9*4* 

12 

IF (P2S 115, 15, 1 ! *4 

00555 

1 95 • 

1 1 5 

IF ( P3 ) 15, 15. 1*4 

00560 

196* 

I 1*4 

IF (P3) !*4, 1 5 « 15 

00563 

J97* 

1 *4 

NLXI ILOC)- M Y 2 

0056*4 

1 9 fl • 


W L Y ( ILOC)- NXI 

00565 

199* 


NLX ( ILOC ♦ 1 )* Nyt ~ 

005 6 6 

200 * 


NLV( ILOC ♦ 1 )- MX 2 

00567 

201 * 


WRITE (KWR, 30) P2‘, P3, NXf, NY | 

005 7 5 

' 02 * 


ILOC- ILOC ♦ 2 

00576 

^03* 


IF (ILOC - IS ID JO, 5, 5 

0060! 

20*4* 

15 

IF (PI ) 110, 10, 116 

0060*4 

205 • 

1 16 

IF (P‘4) 16, 1C, 1 3 

00607 

206* 

1 10 

IF fP*t) 10, 10, 16 

006 1 2 

207* 

! 6 

NLX(ILOC)* nY| 

eck 1 3 

208* 


nlymloo- MXl 

006 1 ** 

209* 


NLXMLOC ♦ i )* NY 2 

00|6 1 5 

2 10 * 


NL Y ( ILOC ♦ f )« M X 2 

006 t 6 

21 I * 


WRITE (KWR, 3D Pi, P *4 , N X |, NY 1 ' 

0062** 

2 l 2 * 


i ILOC « ILOC *2 

00625 

2 t 3 * 


IF ( TLOC-IS 7 ! ) 10 , 5. 5 

00630 

21*4* 

»0 

CONTINUE 

606 3 2 

2 15* 

9 

continue 

006 3*4 

216* 

! ' 

IF (ICT - 8 ) 20, 201 , 7 

66637 

2f7* 

20 

00 119 JJ- 1,16 " 

006*4 2 

218* 


00 1 ?0 K<* l , 128 

006** 5 

219* 


Y( JJ, KK>- y(JJ*) 6 , KK) 
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00606 

220* 


120 

006SO 

221* 


119 

006S2 

222* 



00653 

2?3* 



00659 

2 20* 



00660 

2 25 • 



00661 

226* 



00660 

227* 



00667 

220* 


230 

006 7 3 

229* 



00670 

230* 



0C67S 

23 19 


i 

00700 

2 32* 


j 

00701 

233* 



C0702 

230* 


731 

0070H 

235* 



..... 00705 

236* 



00710 

2 37* 



007 | 1 

238* 


?32 

007 I 3 

239* 


729 

007|6 

200* 


3!93 

00772 

29 I* 


390 

QitZZ.3_„. 

202* 



00726 

203* 


i 

007 3 1 

200* 



00732 

205* 



00733 

206* 


166 

007.35 

207* 


165 

tC2J5 

200* 

C 

i 

00737 

2.09* 



00792 

250* 


i 

007 0 3 

251 • 


i 

00709 

252* 

0 


90706 

253* 


1 

_C02.9 6_ . 

259* 

t 


007 0 7 

255* 


j 

00750 

256* 


?r\ 

0075 3 

257* 


j 

007 5.4 

758* : 



CC757 

259* 


??o 

0076 1 

260* 

719 

00763 

26 1* 


! ' 

00769 

262* 



0C765 

263* 



00766 

?60* 



0076 6 

265* 

C 


007 A 7 

266* 


7 

— eo7fo ' 

267* 



00773 

268* 


13? 

00773 

269* 

c 1 


00779 

270* 



OC 7 7 7 

27l • 



01000 

272* 



Traci 

273* 

1 R n 1 

C 1 0 0 1 

270* 

C 


01003 

275* 




continue 
CONTINUE 
1C • 1C ♦ t 

CAt-L REDTpR<8J 2% I ER i NO, j, IC^T, 70<4fl, TE**P|I*> 

WRITE IKWR, 909 ) NUt 1CNT 
ICT • ICT ♦ T" “ " ' " 

WRITE IKWR, 908) |CT 
IF INW-2009) 230* 229, 230 
WRITE (KWR, 392) NW , IfNT 
KK« 20H9-NW ♦ ) 

t* I . 

00 2.31 K * KK, 2f>09 
TEMPI*). TEmP(J) 

J» J ♦ I 
CONTINUE 
K2» KK-1 

DO 232 *■ I, *2 __ _ 

TEMP IK )• 0. 

CONTINUF 

IF HER -U .393, 399, 393 

WRITE (KWR, 397) ICR, NW 

1 ■ I 

00 t6S JJ» i 7 | 32 

00 186 K*» |, 12« 

V ( J J , KK)» TEMP ( t ) 

!■ 1 ♦ 1 
continue 
continue 

REMOVE STRlpE __ 

00 005.7 IN* 17, 32 

!*■ I IM-I ) • ) 28 ♦ III 
T£MPI IX )- TpMPI I X ♦ I ) 

CONTINUE 
1 1 ■ I 

GO SFaRCH FjRST COLUMNS *GA|n 
GO TO ?l I - -- 

00 219 JJ. ) , 14 

00 270 **■ ) , I2B 
V IJJ, KK)« y I J J+ I 6 , KK ) 

continue 
CONTINUE 
« IC ♦) 

ICT* 10 
IIP » 

GO TO > ; I 

this completes search for «2fro gradient* points 

1 F IN* 1L0C m I 
WRITE (K'lR, (32) J F I N 
FORMAT IIH , 5H t E 1 N» , T ) 0) 

INVERT X ARRAY FOR PLOTTING 
0 0 1 R 0 1 I • I , I F I N 
NX* *'L X C I ) 

MtXlI |» IKAXM-MX 

CONTINUE 

PLOT TAPE 

call GRioivio.')«n,t7e.rt»o.o»i2fi»o» , *»o,R*o,-s,s l -iij,-i(:«*‘‘»* ll *i 
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01009 

276* 


I F I « 1L0C-2 


oions 

27 7* 


DO 800 I« ! , !F| , 2 


OfOlC 

27 8* 


X- NLX(I) 


0101 1 

27 9* 


CALL XSCLVKX, NXp, I ERR ) 


01012 

2 80* 


NX 1 ■ NXP 


01013 

28 1* 


X|« NLX M*1 ) “ " 


0 1 0 1 H 

282* 


CiALL XSCLVKX, MXP, I ERR ) 


01015 

2«3* 


NX 2« NXP 


01016 

289* 


TY* I NMAV+I )«NLY< I > 


01017 

285* 


CALL YSCLVKYY, NyP, IfrRj 


01020 

286* 


NY I ■ NYP 


”5)1 C 2 i 

'287* 


YY* CNMAY+i )-NLV< 1*1 1 " 



01022 

1 

288* 


CALL YSCLVKYY, N Y p # IFRR) 


011023 

2«9* 


N Y 2» NYP 


Ojt 02 M 

29Q* 


CALL LINFVINXI, NYI, NX2i My2> 


oil 025 

29 | * 

803 

CONTINUE 


0,10 27 

292* 

999 

continue 


01 03 1 

293* 


CALL ENDJOR ' 


0103? 

2 99* 


STOP 


01033 

295* 


END 



TNO OF 

COMPILATION? no DIAGNOSTICS. 



^ I ^0^ J J‘ AG E IS 

°° R fiVA2JXr i 


FIGURE 2 - Concluded 
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MAGNETIC FIELD GRADIENT PROGRAM 

One way to represent magnetic field energy distributions in 
sunspots is to plot values of magnetic field gradient as a function of 
position in the spot. A second special-purpose program has been 
written to generate such displays. 

The value actually plotted is not the value of magnetic field 
gradient; instead, it is the difference between signal values of adjacent 
grid positions. This procedure was chosen for several reasons. First, 
it | gives the most direct representation of the actual data. Second, the 
plotted values are independent of uncertainties in the field strength 
calibration of the signal. Third, evaluation of distances on the solar 
surface and foreshortening corrections are not required, so the plot 
is independent of the assumptions which must be introduced to make 
those calculations. Fourth, if differences only are used, the evaluation 
of a numerical derivative is avoided; that procedure is notoriously 
unstable. The final plot is presented in the original magnetogram grid 
format, so it can be easily compared with contour plots and with photo- 
graphs of the scope display at the time of observation. A label page 
giving the date and time of the observation is written for each plot. 

Difference values in three size ranges are plotted. The scale 
factor used for demonstration was the approximate mean noise level of 
the observations, XNSE, but any interesting scale factor could be 
substituted. The plotting symbols were chosen to give an appearance 
of increased density as the signal difference value increased. The 
difference ranges and corresponding symbols are listed in Table 1. 

The magnetograph signal data are read from magnetic tapes in 
the same format as are used for the neutral line program. In addition 
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TABLE 1 . GRADIENT PLOT PROGRAM SYMBOLS 


GRADIENT DIFFERENCE 

SYMBOL 

| 2 • XNSE £ difference < 5 • XNSE 

0 

5 • XNSE 5 difference < 10 • XNSE 

+ 

difference 2 . 10 * XNSE 

* 


to the data tape, a mean noise level or other plotting scale factor must 
be specified. This value is supplied in a DATA statement by assigning 
a value to a variable called XNS; the appropriate noise value, XNSE, is 
later calculated from XNS. In the sample given here, XNS = 0. 05. As 
in the neutral line program, the number of matrices to be processed is 
supplied in a DATA statement. Again, the variable name assigned to 
this quantity is IPP, and again it is an integer variable. This program 
processes either longitudinal or transverse data but cannot do both types 
in a single run. If longitudinal data are to be plotted, the magnetograph 
data tape, after initial processing, is used as the data source. IPP is 
the number of longitudinal matrices on the tape. A third DATA state- 
ment must be included: 

DATA LT/1 / 

and will cause the program to select the correct procedure for reading 
and plotting the longitudinal data. If transverse data are to be plotted, 
the data tape must be generated by the fourth program described in 
this section, IPP is the number of transverse matrices to be plotted, 
and LT must be supplied by the DATA statement 

DATA LT/2/ 

No separate data deck is required. 

The logic of the program is like that of the neutral line plot 
program, except that value differences instead of sign differences are 
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recorded. Sequences of small grids, numbered PI, P2, P3, P4 as 
shown in Figure 1, are used to evaluate the differences across the 
magnetogram. The plot symbol is assigned to the coordinates of either 
PI or P2. The coordinates of PI are used if one of the differences 
(P1-P2), (P1-P3), or (P1-P4) is plotted; the coordinates of P2 are 
used if (P2-P3) or (P2-P4) is plotted. The value plotted is the largest 
difference found in the small four-point grid. No difference values 
smaller than 2*XNSE and no isolated difference values larger than 
6*XNSE are plotted. 

A listing of the program is given in Figure 3, and a sample 
plot is shown in Figure 5. 
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00100 

oomo 

00100 

ootoi 

00103 


001 OS 
00106 
00107 
001 10 
001 I I 


BLANK J S PLOTTED IF GRADIENT VALUE IS LESS THAN 2X NOISES 

o if gradient is 2 to sx noise? • if gradient is more than 
10X noise. ALL GRADIENTS MORE than 2X noise ARE PRINTED. 
OlMEMSiflN X TWO < I 000 ) , VTWO ( 1000) , Xp l WE ( I 000 ) . VF I VE < 1000) 
DIMENSION XTEN(SOO) lYTEN(SOO) .TEMPI5000) ,Y<32, | 2 B ) 


DIMENSION FLOX ( 1 2 ) , FLDYI12. 
DIMENSION A ( 1 9 ) 

0 I MENS I ON L T I ME ( 9 ) 

DIMENSION TEME 12099) 

DAYA IPP/38/ 


00 M 5 
00117 
00121 
00123 
001 25 
- & 6 4 - 2 S - 

00127 

oom 

00132 

00133 
00139 


00136 

00136 

00137 
001*0 
00191 


DATA FloX(2)/6H / 

DATA (FLDXUli 1« 9, I2)/9«6H / 

DATA (ARRAY(I>, I ■ l , 22 ) / 2 9 H HARDCOPY ONLY, ONE COPYc 18«6H 
DATA HL/8/ 

data LT/1/ 

IT . ) I P L ONGI T U DI N AL . 3 I F T -BA-NS- M fcRS-C . 

DAJA XNS/.05/ 

KHE ■ S 
K«R» 6 

CALL OPEN ( 8 * 1, S) 

IMAXb 128 


1S1ZE « 1000 

NO, OF ROMS, NO. OF COLUMNS, NO. OF GRADIENT VALUES STORED 
■151b loco 
1 S2» 1000 
IS3b SOO 


00170 

00170 

00171 

00172 

00173 


CALL 1DENT (9, ARRAY) 

DO 99 7 I P b 1, |PP 

I read header information 

300 CALL KEdTPR (8, 2, | £R , W* r 19, A) 

j IF ( N A » 1 9 ) 300, I, 300 

j990 FORMAT ( IH , 3 HNWb, IS ) 

IF I IP-30) 997, 700, 700 
700 CONTINUE 

IF ( L T- 1 ) 1 70, 1 70, 171 
;170 X N 5 E b A { 7 ) • ( XNS- A ( 13)) 
j -fro- - T O— 1 7 2 

171 XNSE* SORT I A ( B ) *S(JRT ( XNS ) ) 

XNSE IS MEAN NOISE LEVEL 

172 TXNS> XNSE.2.00 
FXNSb XNSE*5 .00 
DXNSb XnSE* 10.00 


FIGURE 3. LISTING OF MAGNETIC FIELD GRADIENT PLOT PROGRAM 
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00201 
00205 
0020i 
002 t H 
0021 5 


391 WRITE (KWR, 392) IER. NW 

392 FORMAT t I H , 21 1 0 > 

390 WRITE (KWR, 4.7) <A<I), I* 1, 19) 
67 FORMAT IIH , 61 1 X . 012) ) 

WRITE (K*R, 167 ) 1 A< 1 ) . 1 • 1 , R ) 


0022R 
0023 2 
00233 
002R1 
002,R 1 

58# 

59# 

60# 

61# 

62# 

169 

1 6B 
C 

WRITE (KWR, 169) <A(I), t* 5, 11) 

FORMAT <1H , 5 ( E 1 2 , 5 , 2 X ) , 2 ( 1 1 0 , 2 X ) ) 
WRITE (KWR, 168) ( A C t ) , 1* 12, 19) 

format uh , e(Ei2.5, 2x>) 

READ 1ST 2 RECOROS OF MA&NETOI.RAM 

0o2iR5 

6 3# 
6R* 


l/U I ■ 1 | i ' ... .. 

ICT * I 

0O2R6 

65* 


WRITE (KWR, 9881 ICT 

00 25 1 

66# 

988 

FORMAT < 1H .RHICTc, 15) 

00252 

67# 


J* 20R8# ( l-l ) * 1 

00253 

68* 

19 

CALL KEDTPR <8, 2, IER, Nw, 1, ICnT, 20R®, TEMP(J)) 



00276 
00301 
0030? 
00303 
00305 
-MHO 
00317 
00315. 
00315 
00317 
0032? 
003 2 3 
0032R 
00326 
00331 
00335 
003R0 
— 0 03 R) 
003H? 
003R3 
003MR 
003 RR 
003RR 
-06 3 4 5 - 
003H6 
003 R7 
003R7 

00350 

00351 


989 FORMAT ( I H , 3HNW«, 17, 5X, 5MICNT*, 15) 
IF (NW-20R9* 130. 129, 130 

130 WRITE (KWR, 392) NW, ICNT 
00 122 | B- 1, 3 
122 TEME (IB)* TEMPI IB) 


00 123 IB* 11, 20R9 
TEME IIBI> TEMP(IMB) 

I 1MB* 1MB ♦ t 
123; CONTINUE 

DO 1 2 R JB- R , 10 

00 125 IB* 1, 20R9 
1 25 TEMP( IB ) • TEMEC IB ) 

REMOVE STRIPE 
129 00 R056 I M« 1 , 16 

IX* ( 1M.1)« | 28. J-l ♦ 1 1 1 

R056 CONTINUE 

IF (IER-1) 3’5, R650 , 395 
395 WRITE (KWR, 392) IER, NW 
R650 IF ( 1-1 ) 65, R 65 , 65 
R65 lx • ion 

1 DEL V a 25 

N* FLU( 7, 10, TEMP( 1) ) 

10 AYS* 200 • F Lo (26,1, N)*I00*FLD(2?,1,N)»80*F LO <28.1, Nl 
1 ♦R0*FLD ( 29 , 1 ,N)«20»FLD ( 30 , l ,N 1 + 1 0#FLD ( 3 1 . I ,N ) #a#FLD ( 32 , 1 ,N 1 
2*R*FL0(33,l,N)+2*FLD(3R,l,N)*FLDl35,l,N) 


FL0(30,R.N)» FLD( 12,R,TEMP( 1 1 ) 

1 HOURS* 20#FLp(3Q l l,N)*t0#FL0(3l.l.N)#8#FLD(32,|,N) 
1>R*FL0( 33, 1 ,N) + 2*FLD( 3R , 1 ,N)*FLO <35, 1 ,N) 

N* FLO( 20, 7, TEMPI ] ) ) 

1 M I N* R0*FLD<29,1 , N ) ♦20*FLD < 30 , l , N ) ♦ 1 0#F»-D ( 31 . 1 , N ) 


FIGURE 3 - Continued 

ORIGINAL PAGE IS 
OP POOR QUALITY 
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003*2 106* 

00 ** 3 1 0 7 * 

003*3 108* 

00354 109* 

00362 ltO* 

00363 ill* 

00369 112* 

00*6 3 f l- 3 * 

00366 119* 

00370 115* 

00370 116* 

00370 1 1 7* 

00370 na* 

80*70 h- 9 -*- 

00371 120* 

00372 12|* 

00376 122* 

00400 123* 

OOfOl 129* 

00609 126* 

00909 127* 

00006 l 2 fl • 

00907 129* 

00910 130* 

009| 1 m*- 

00912 132* 

00912 133* 

00j9 I 3 139* 

00913 135* 

00919 136* 

0 04 1 7 - +37*- 

00927 138* 

00923 139* 

00926 190* 

00927 191* 

00932 192* 

004 +J - -+ 40 * 


00939 
00936 
00990 
0099 1 
00999 

nn ii a i 

144* 

145* 

146* 

147* 

148* 

uu '“o 

00950 

1 4 9*" — - 
ISO* 

00950 

151* 

00950 

152* 

00950 

153* 

00950 

154* 

--0046* 

-155*- — 

00964 

166* 

00965 

157* 

00965 

158* 

00960 

169* 

00961 

160* 

-004+2 - 

Hr**- 


N* FLO ( 27 , 7 , TEMP ( 1 ) 1 

U E C * -4 Q « FL - p l3 * .l .N)*-20* F L&<30,1 , 14 > * 1 0* FL P ( 3 1 >1 , N > * fl >F LP t 32 . 1 ,N) 

1*9*FL0133,I,N)*2*FLD(39,1*N)*FLD«35,I.N) 

ENCODE (90, LTIME) lOAYS, l HOURS ♦ ININ, I SEC 
90 FORMAT (19, SH DAYS, 17, 6H HOURS, 18, 8H MINUTES# |8,8H SECONDS) 
CALL FRAMEVO) 

CALL PR1NTV (59, LTIME, IX, 1 V > 

CA L L STp PT V — 

65 CONTINUE 
I C« l 

C 1 C T IS NO, of RECORDS READ FROM TAPE 

C 1C 15 NO, OF TIMES READ FROH TAPE IS EXECUTED 

C 1C should be equal TO l C T - | 

-€ O ftflE - ts T -g- SECO B B S + N TO - Y -6 R RA Y 

I* 129 

00 68 JJ* 2, 32 

DO 70 KK* 1 , 128 

YC JJ. KK I- TEMPI I I 

1 ■ I ♦ 1 

68 CONTINUE 

C SET X COORDINATE COUNTER AND START SEARCH OF FRAME 

NM A X » 1 MAX -I 
Jl* 1 
J 2 * 1 


U3*-l - - 

11*2 

C START SEARCHING COLUMNS 2 THROUGH |6 

NM A Y * IMAY - l 
C SPIKE NOISE FILTER 

2 1 1 DO 9050 JJ* I , 31 

P0 - 9Q S1 K K* 1 , I 2 7 


XTEM* ABS(XNSE-A05(Y(JJ,KK) I ) 

IF (XTEM-3.*XNSE > 9051, 9051, 9052 
9052 Y T E M ■ ABS(XNSE-ABS( Y< JJ+l ,KK) )) 

IF ( YTEM-3,*XNSE ) 9053, 9053, 9051 
9059 Y ( JJ ,KK 1 * XNSE 
fr O — TO " 405-1 — 


9053 2TEM* A B S ( XNSE- A 9 5 ( Y < J J , KK ♦ 1) I) 

IF (ZTEM-3.*XNSEI 9065, 9055, 4(}S| 
9055 ZZTEM* A9S!XNSE-AeSlY(JJ*i,KK+lll) 
IF (ZZTeM“3o*?NSE) 9059 , 9059, 9051 
4061 CONTINUE 


00 j 9 JJ* 11,16 

PREPARE TO GET CORRECT COORDINATES FOR (jRAQIENT VALUES TO 
8E STORED 

MISUNDERSTANDING about ORDER in H»h!CH tape AAS written 
SO MUST EXCHANGE X AND Y COORD AFTER THEV ARE DETERMINED 

*~ <4-C« | >*» 6 - 

NX 2* N X 1 ♦ 1 
DO 10 I- I, 127 
SET Y COORDINATE COUNTER 
NY l * 1 
NY 2* 1 ♦ I 


FIGURE 3 - Continued 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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00463 

QQ444 

16 2* 

1 

P2» Y ( Jj, I ♦ 1 ) 

P2* Yt-JJ-* 4-, l_) 

00465 

164* 

P4» Y ( J j* J , 1*1) 

00466 

165* 

GTEHU )■ PI - P 2 

00467 

166* 

GTEM ( 2 ) ■ P1-P3 

00470 

167* 

GTEM ( 3 > «P 1 *»P4 

00471 

— ww — 

168 * 

GTEM ( 4 ) » P2-P3 

Hr9* 

GT-E4* (-5 — P*»P-4 


J 70* 
I 7 1 • 
172* 
173* 
174* 
1- 7 5* 


00525 
005 26 
00527 

00530 

00531 
00&3 4 


182* 
18 3* 
|88* 
185* 
l 86* 
—4-8-7-*- 


00537 

00544 

00545 

00546 
00557 

- 0 05 57- 
00557 

00560 

00561 
0056? 
00565 


00570 

00573 

00600 

0060! 

0060? 

-00^3- 


188 * 

189* 

1 90* 

1 9 J • 

» 9 2 • 
4-9**- 
194* 
195* 
196* 
197* 
198* 
- > 99* 
200 * 
201 * 
202 * 
20 3* 
204* 


00604 
00607 
006 ! ? 
006j|7 
00620 
-Oft *? -+ - 
00622 
00625 
006 32 

00633 

00634 


2 06* 
207* 
208* 
209* 
210 * 
-2-H-*- 
212 * 
213* 
214* 
2 15* 
216 * 


301 


GMAX« A0S<GTEM( 1 I ) 

IG • 1 

DO 299 12- 2 i 5 

IF <ABS<GTEM< 1 2 ) l-GMAX) 299, 301 * 301 
GMAX* GTEM ( I 2 ) 

-TG* A i 


00473 

00474 

00475 
00500 

00503 

0 05 0 4 

00505 176* 299 CONTINUE 

00507 177* IF (16-4) 297, 29», 298 

00512 178* 297 IF (GMAX-TXNS) 1 0 , 295, 295 

00512 179* C COORD OF PI WILL BE STORED 

00515 180* 295 IF <GMAx-FXnS> 294, 302, 302 

- 00 5* 8 4-8-4* 2 8 4 W 8 H E K WR , 2 % 4-- N -H A U l , GMAX 


25 FORMAT (IH , 4HN X 1 ■ , I 7 , 9H 
XTWO < J ] ) « Ny 1 
V T WO < J 1 >* N X 1 
J 1 » J 1 ♦ 1 

IF (Jl-ISM 10, 5, 5 


NY 1 * 1 1 7 , 9H 


GNAX-,E12.5) 


303 WRITE IKWR, 251 NX1, N Y 1 , GMAX 
XFIVE1J2)* NYl 
YF ! VE I J2 » * NX 1 
J2» J2 ♦ 1 

j IF (J2-JS2) 10, 5, 5 

. .biXL»_GMAJ( 

! X T E N < J 3 ) • NYl 

YTEN|J3»» N X 1 
J 3 » J 3 ♦ 1 

j IF <J3-IS3> 10, 5, 5 

298 IF (GMAX-TXNS) 10, 307 , 307 

-COOWD—OF-P2- W 1 LL 

307 IF JGMAX-FXNS) 308, 309, 309 

308 WRITE (KWR, 26) NX2, NYl, GMAX 

126 FORMAT ( I H , 4HNX2« , I 7 , 9H NYI»,|7,9H 

XTWO ( J I I* NYl 
YTWO(Jl)* NX 2 

-!■ 4 -- 


GMAX ■ • E 1 2 • 5 ) 


309 
3 10 


311 


IF < J 1 - I S | ) 10, 5, 5 

IF (GMAX - DXNS) 310, 311, 311 
WRITE (KWR, 26) NX 2 , NYl, GMAX 
X Fi 1 V E < J ? 1 * NY 1 
YF,1VE<J?>« N X 2 

- — (J2 ■*• 4 — — — — > - — - - — - 

IF (J2-IS2) 10* 5, 5 

WRITE (KWR, 26) Nx 2 » NYl, GMAX 

X T E N < J 3 ) * NYl 

Yt£N < J3 ) * NX 2 

J3> JJ * 1 

-4#— 4**~4$3> 10, 5, 


ORIGINAL 
.015 POOR QUALITY* 


FIGURE 3 - Continued 
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00**40 

0 0**4 3 - 


006 *4 6 
00**47 
0065 *4 
006 5 i 
00660 
0 066 1 


00662 

00662 

00665 

00666 
00672 
00 6 7 3 
00676 

006 7 6 

007 0 I 

00*0 3 

00705 

0fr*- h 0~ 

0071 3 

007 16 

00717 

00721 

00723 

007 2 -* 4 - 

007 2*4 

0Q725 

00730 

00731 

00732 

nn7 

Uw r 

007*40 

007*4 1 

007*4*4 

007*47 

00753 

©09**4- 

00755 

00760 

00761 

00762 
0076*4 


218* 

21 9* 

220 * 

2 21*- 
222 * 
223* 
229* 
2- 2 5*- 
226 * 

2 27* 
228* 
229* 
230* 
231 * 
232* 
233* 

2 3*4* 
235* 
236* 
237 *- 


5 IF I NX | - 70 > 6, 6. 906 

6 W RITE (K »H 1 * ) * * 4- 


l 1 

906 

907 


FORMAT <6m NYI«, 17) 

WRITE fKWR, 27) Jl, J2, J3 
WRITE (KWR, 11) 

WfMTE UtWR , 12 ) 

format c*«oh array reserved for gradients is filled.) 

FO RM AT ( 10 H- TE R-H l N»T€ -ftWr ) 


IF <ICT - 8) 907. 7, 7 
GO TO EN& OF MAGNETOGRAM 

CALL REDTPR (8. 2, JER, NW, 1. ICnT. 20*4». TEMP { 1 ) ) 
WRITE +KWR, 9-89 ) NW , ICNT 
ICT ■ ICT ♦ 1 

- # R1 TE 4K WH . 9 88 ) ICT 


IF IICT-8) 907, 7, 7 

GO TO END OF MAGNETOGRAM, THEN START NEXT ONE 
10 jCONTlNUE 
9 CONTINUE 

IF (ICT - 8 ) 20, 201 , 7 
■ 20 DO 11 9 ) , — 1-6- 


007 7 n 
00771 
007 7 3 
00776 
01002 
0 1 0 0 3 ' 
01006 
01011 
01012 
01013 
01015 
oro+9 


23e* 
239* 
2*40* 
2*4 1 * 
2*42* 
-2-*4- > *~ 


2*4*4* 
2*45* 
2*46* 
2*47* 
2*48* 
2 * 49* 
250* 
25J* 
252* 
253* 
25*4 • 
2fr&*- 
256* 
2 57* 
25 8* 
259* 
260* 
■ X 6 f *~ 
262* 
263* 
26*4* 
*65 • 
266* 
2 6 7 * 


00 120 KK« 1 , 128 
Y( JJ, K K ) ■ Y( JJ ♦ 16, KK) 

120 CONTINUE 
1 19 CONTINUE 

1C - 1C ♦ 1 

C A LC. — R E 0 T P-R — ( 8 , — 2 N#, — 1 r — T-CN T , -20*48, T E M P 4 -I.)) — 

REMOVE STRIPE 
DO *4057 IM* 1 , 16 

IX* ( 1M-1 ) • 1 28* 1 1 1 
TEMPI 1 X )« XNSE 
*|C57 CONTINUE 

ICT ■ ICT ♦ l 
WRITE (KWR. 988) ICT 
IF ( N 6-20*49 ) 230 . 229, 23 q 
230 WRITE (KWR, 392) NW, ICNT 
KK* 20*49 - NW ♦ 1 


231 


00 231 K ■ KK. 20*49 
TEMP ( K > « TEMP(J) 

J* J ♦ 1 
CONT 1 NUE 
K2» KK - 1 


268* 
269* 
2 70* 
271* 
272* 
27> »- 


TEMP ( K ) a 0, 

232 CONTINUE 

229 IF (ItR- 1 ) 393 , 39 S, 393 
393 WRITE (KWR, 392 ) IER, NW 
399 I« 1 

DO 166 kK> l , 128 
Y( JJ, KK)» TEMPI I ) 

I ■ I ♦ 1 
166 CONTINUE 
1 65 CONTINUE 

j-H* — f - • — • -- 


FIGURE 3 - Continued 
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RESULTS OF PLOT PROGRAMS 

Figures 4 and 5 show results of the two plot programs for the 
same frame of data (23 October 1973, 23^ 41 m 9 s U. T. ). Figure 4 
is, the neutral line plot and Figure 5 the field gradient plot. Figure 6 
is| a contour plot of the same data frame, made by the Computer 
Sciences Corporation. (The orientation of Figure 6 is reversed both 
vertically and horizontally from Figures 4 and 5. ) The contour 
intervals of Figure 6 are labeled according to the code given in Table 2 

TABLE 2. CONTOUR LABELS FOR FIGURE 6 i 


I The position of the solar limb shows clearly in each of these 
three figures. (The vertical streak in Figures 5 and 6 is a flaw on 
the magnetograph vidicon tube. The gradient plotting program now 
contains a routine which removes this streak. ) In addition, the active 
region shown contains two neutral lines which run across the regions 
of steepest magnetic field gradient. There is also some evidence, 
in Figures 4 and 6, of a weak but complex field structure near the limb 


LABEL 


FIELD STRENGTH (gauss) 


A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 


-21581 

-1.290 

-6.452 

-2.581 

-21581 

0.000 

+2.581 

+2.581 

+6.452 

+11290 

+2.581 


x 10 3 
x 10 3 
x 10 2 
x 10 2 
x 10 1 

x, 10 1 

x 10 2 
x 10 2 
x 10 3 

X 10 3 
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FIGURE 4. SAMPLE NEUTRAL LINE PLOT 
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FIGURE 5. SAMPLE GRADIENT PLOT 
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FIGURE 6. MAGNETIC FIELD CONTOURS FOR DATA PLOTTED 
IN FIGURES 4 AND 5 
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NOISE REMOVAL PROGRAM 


Some of the data matrices contain noise values which appear 
as single values much larger than any other values near them. This 
program removes each of these values, replaces it by the average of 
the points surrounding the removed value, and writes a new tape 
containing the corrected matrices in the same format as on the original 
tape. Both longitudinal and transverse matrices can be corrected in 
one pass. 


The program also removes the vertical streak which appears 
in Figures 5 and 6. This stripe always occurs in column 1 1 1 of the 
data matri: as it is read in (the matrix is stored reading from right 
to left across the frame). Because the location is always the same, 
the stripe can be removed as soon as the matrix has been read. The 
value in the location X = 111 in each line is replaced by the value 
X = 112 in the same line. 


The isolated values to be removed are identified by comparing 
them to criterion values which are supplied as data. Two criterion 
values are needed, one for longitudinal data and one for transverse. 
The data tapes used are the magnetograph data tapes after the initial 
processing, so the data values are signal values, not field strengths. 
Also, the transverse criterion value is the criterion for the U or the 
R component, not for the resultant transverse field vector. It is 
assumed that the same criterion value can be used for the U and R 
matrices. 

After a matrix is read, the proper criterion value is chosen, 
and the absolute value of each matrix data point is compared to the 
criterion value. If the matrix value is more than six times the 
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criterion value, the adjacent matrix points are examined. If one or 
more of the adjacent values is greater in absolute value than six times 
the criterion, the tested value is assumed to be valid. If none of the 
adjalcent values is greater in absolute value than six times the criterion 
the tested value is assumed to be invalid and is replaced by the average 
of the adjacent values. 

The magnetograph data for the program are taken from the 
data tapes after initial processing. These tapes have the general 
format shown in Figure 7. 


MAGNETOGRAM 

#1 


HEADER 

[LONGITUDINAL DATA MATRIX 
INTENSITY DATA MATRIX 

U MATRIX (ONE COMPONENT OF TRANSVERSE DATA) 

R MATRIX (OTHER COMPONENT OF TRANSVERSE DATA) 


i HEADER 

LONGITUDINAL DATA MATRIX 
INTENSITY DATA MATRIX 
U MATRIX 
R MATRIX 


(y 


HEADER 

J LONGITUDINAL DATA MATRIX 
MAGNETOGRAM < INTENSITY DATA MATRIX 
#NMAT / U MATRIX 
R MATRIX 


V F 


END OF FILE 


FIGURE 7. MAGNETOGRAPH DATA TAPE FORMAT 
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On any tape, the longitudinal matrix may or may not be present, 
the intensity matrix may or may not be present, and the U-R matrix 
pair may or may not be present. This program will process any 
combination of matrices if the correct data are supplied. (It is 
assumed, however, that the intensity matrices are either always 
present or always absent. All tapes processed so far have complied 
with this assumption. ) 


In addition to the tape, values must be furnished to the program 
through four DATA statements. The contents and formats of these 
statements are described in Table 3. 


TABLE 3. DATA STATEMENTS FOR SPIKE 
NOISE REMOVAL PROGRAM 


i 

VARIABLE 
j NAME 

VARIABLE DEFINITION 

MODE 

(REAL OR 
INTEGER) 

EXAMPLE 
OF DATA 
STATEMENT 

NMAT 

Number of magnetograms 
on tape (i.e., number 
of groups of matrices 
as shown in Figure 7) 

Integer 

DATA NMAT/ 7/ 

XNA 

j 

Criterion value for 
longitudinal data = 
six times largest 
valid isolated signal 
value ! 

. . j i 

Real 

DATA XNA/. 005/ 

XNB 

i ■ 

i. ... ; ' 

i r ■ ■ 

i ; 

: 

i, _ i' 

Criterion value for 
each component of 
transverse data = 
six times largest 
valid isolated signal 
value expected in U 
and R matrices 

Real 

■ ■ ■ 

DATA XNB/. 002/ 

DINT 

| 

i 

DINT - 1.0 if intensity 
matrices are present 
on tape and =0.0 if 
no intensity matrices 
were recorded 

Real 

... 

DATA DINT/0.0/ 
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After the matrix has been corrected, the program writes it on 
tape in the same format in which it was read and reads the next matrix. 
The tape which is produced contains NMAT groups of matrices; each 
group contains the same number of matrices as were in that group on 
the original tape. 


A listing of the program is given in Figure 8. 


00 100 

1* 

C 

KEMOVE spike NOiSe. values 

00101 

2* 


01 MENS ION TEMPI 16 38H 1 , A < 1 5 > , I A I H 1 , Y ( 1 28 . 1 28 > 

00103 

3* 


DIMENSION TEMEf*C-H9) 

00 1 OH 

H* 


01 MENS.I0JN . XN5.E1 S ) . . 

00 1 OH 

5* 

c 

nmat« no. of magnetograms on data tape 

00105 

6* 


Data nmat/7/ 

0 C l 0 7 

7* 


data xna/. 00583 / 

03111 

8* 

i 

DATA XNB/.Q015A/ 

001 1 1 

9* 

c 

XNA IS LIMIT FOR L MATRIX 

001 1 1 

10* 

c 

XNB IS ^.lALlI..XOa T MATRIX 

00113 

11* 


DATA 0INT/1.G/ 


FIGURE 8. LISTING OF NOISE REMOVAL PROGRAM 


ORIGINAL PAGE Ifc 
OP POOR QUALITY 
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aoi »3 

12* 

C 


DINT-0 IP. NJLLiaLfclSJ TY._Ha.TR 1CLS * *1 If .INTENSITY H*TR If E5 

OOI 13 

13* 

C 


present 

00 115 

1*4* 



KrE* 5 

ociiA 

15* 



KAR* 6 

0C1 17 

16* 



CALL OPEN <8»1,3) 

00120 

17* 



call remind <s> 

00121 

IB* . 



call sj. ... ...... 

00122 

19* 



call rewino<9> 

00123 

2U* 



Do 1 00 IMaT- It nmat 

00123 

21* 

c 


READ HEADER 

00126 

2 2* 


1 

call reotpr < a , 2 » ier ,nw , 2 , i a< n * 7 . a < n . 2 , i a ( 3> , 8 , a<«> > 

00127 

23* 



CALL WRITER <9,2. IER.2, I A < 1 ) ,7,*m ,2, IA53I,8 ,a(8H 

00130 

2*4* 


. ... ... 

.lx JN^-4 91 L». 1.7. . 1 . ,■ 

0C1 33 

25* 


17 

WRITE (KA'R, 21) NW 

00136 

26* 


21 

FORMAT (<4h NW«, 15) 

00137 

27* 


26 

WRITE (KWR, 167) 1 A ( 1 ) , 1 A 4 2 ) • A f 1 ) . A 1 2 > 

001 **s 

28* 


167 

FORMAT (1H *2 ( I 10. IX ) ,2 ( A12.2X H 

001*46 

29* 



WRITE (KWH, 169) (<A<I).1= 3. 7 ) , 1 A < 3 ) , I A ( *4 ) ) 

00 15 6 

3 0 * 


4 69 

Format u h »5iei2.5.2x) » 2 < jio. 2 x) ) 

3C 167 

3 1* 



WRITE CKWN, 16 8 9 ( A 4 I ) , 1 * 8, Ifil 

00165 

3 2* 


168 

FORMAT 4 1 H . 8<El2.5i2X>) 

00166 

33* 



WRITE IkWr, 299f I A ( 2 ) , IMAT 

00172 

3*4* 


299 

format < in . P< I to. sxi > 

00173 

35* 



IF 40INT) 2. 2. 3 

00176 

36* 


2 

IF < I A ( 2.1 -51 *4, 5, 5 

00201 

37* 


*4 

I RUN- 1 

OC-2 J2 

38* 



XNSEH)* XNA 

00 2q 3 

39* 



GO TO 50 

0020*4 

*4a* 


5 

IF < l A ( 2 1 *6 | >6,6,7 

C02o7 

*41* 


6 

IRUN* 2 

30210 

*42* 



Xh|$E( J 1 !• XNB .... 

U02 | 1 

*43* 



X NSE < 2 ) * XNB 

00212 

*4 *« • 



GO TO 50 

0C2 1 3 

*45* 


7 

iron- 3 

0021*4 

*46* 



X N S E ( 1 ) a XNA 

002 | 5 

*47* 



Xr!5E<2>* XNB 

00216 

*48* 



XnSE ( 3 ) « XNB .... ... 

00217 

*49* 



GO TO 50 

0C220 

5C* 


3 

If ( I A < 2 ) - 5) 8. 8, 9 

002 2 3 

Si* 


8 

I RUN* 2 i 

0022*4 

S2* 



XNSE { I ) m xN A 

00226 

S3* 



X NSE 4 2 ) * XNA 

00226 

5*4* 



60 TO 50 - - . . 

0Q|2 2 7 

55* 


9 

if < i a < 2 > -6 1 1 i n • to. 1 1 

0023 2 

56* 


10 

IRUN* 3 

C0233 

57* 



X N 5 P < 1 ) * XNA 

0023*4 

58* 



XnSE <21* XNB 

00235 

59* 



X N S E < 3 ) * XNB 

0023 6 

6o* 



GO TO 50 ... ... ...... 

00237 

6 I • 


1 1 

IRON * *4 

00;2*4 0 

62* 



XNSEH )■ XNA 

0G;2 *4 1 

63* 



*NSE<2)« XNA 

002*42 

6*4* 



X NSE < 3 ) * XNB 

002*4 3 

65* 



X NSE < *4 1 ■ XNB 

002*4*4 

66* 



Gil TO 50 ... ... .. .... 

002*45 

67* 


50 

Do 51 I* 1 , IRON 
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0024 5 
□0250 

. 6 . 8 * — - C- 

69* 

R£ A D MATRIX - — .......... ...... ... . 

J* 1 

002 5 1 

70* 

52 

Call RE0TpR(8,2)IER,NW,I, ICNT ,204 8 , TEMP ( J ) > 

00252 

71* 


write ikwr. 9 89i nw, ic n t 

0026 6 

72* 

989 

Format tin » 3 hnw*. i7,5X,5h!Cnt*, isi 

00256 

00257 

73* C 

Turn 

check for wrong record length in ist record 

—I F- W 2 0 ^ J ^ • 1 0 

■ 

*'“*"** 


0026 2 

75* 

130 

write (kwr, ’89i nw, ic n t 

00266 

76* 


WRITE (KWH t 1130) ( TEMP 1 I L ) » JL« 1* 20) 

□027 4 

77* 

1 130 

FORMAT ( 5 ( 1 H , E12.5)) 

00275 

78* 


00 122 IB* It 3 

00300 

00302 

79* 

Bn* 

122 

TEMEl 1 B) * TEMp ( |B) 

4.MB* 4 - . ... . 



00303 

6 1 * 


DO 123 IB* lit 2048 

00306 

82 * 


TEME(IB)* TEMP(iMB) 

Q03O7 

83* 


IMB* 1MB ♦ 1 

00310 

84* 

123 

CONT INUE 

003 1 2 

86 * 


DO 124 IB* 4 t 10 

00315 

86 * 

1-2 4 

TE-MF( 1B>* XNSE141 — - 

003 l 7 

87* 


DO l 25 IB* It 2048 

00322 

88 * 

125 

TEMP( IB)* TEME< IB) 

3032* 

89* 


WRITE (KWR, 113*) ( TEMP ( I L ) t IL* It 20> 

00332 

9 o* 

59 

J* J ♦ 2048 

00333 

9i* 

53 

Call REDTpR ( 8 , 2 . 1 ER ,NW , 1 , 1 CNT .2048, TEMPI J) ) 

00334 

92* 


WJUTE (K 9.89).. NW* -IcNT 

00340 

93* 


IF (ICNT- 8 ) 12. 129, 12 9 

0034 3 

9 4* 

12 

J* J ♦ 2048 

0034 4 

95* 


60 TO 53 

00344 

96* 

C 

REMOVE STRIPE 

00345 

97* 

129 

DO 4056 I M= 1. J 2 8 

00350 

9 8 * 


IX* ( 1M-1 )• 128*4 1 1 

0035 1 

99* 


TEMPI IX )* TEMP 1 IX ♦ 1 ) 

00352 

100 * 

4056 

CONT I NUE 

00352 

1 0 i • 

C 

WRITE MATRIX INTO Y ARRAY 

00354 

102 * 


I Y* 1 

00355 

103* 


00 68 JJ» 1, 128 

00360 

104* 


DO 70 KK* It 128 

0036 3 

1 Os* 


Y< JJ.KK)* TEMPI I Y ) 

00364 

106* 


I Y * I Y ♦ I 

00365 

107* 

70 

Continue 

00367 

108* 

68 

continue 

00367 

109* 

C 

Spike noise filter 

0037 1 

1 10 * 

211 

Do 4050 JJ= 1, 1.2 7. 

00374 

111 * 


Di) 4 05 1 KK* It 1 2 7 

C>C377 

1 1 2 * 


XTEM* ABS(XNSE( I )-ABS( Y( JJ,KK) ) ) 

oc^cto 

113* 


X MG* Y ( J J , K K ) 

TH>4oi 

1 14* 


IF (XTEM-5»*XNSEU ) ) 4 o5 1 , 4c5l, 4052 

00404 

i i 5 • 

4 05 2 

Y T EM* ABS(XNSE(!>-ABS(YUJ+1,KK>)> 

oc4os 

116* 


IF tYTEM-b,.*XNS6<4 > ) 40&3, 40&3, 4051 

004 1 0 

117* 

4054 

ttmi* ( Y ( J J . k K ♦ 1 ) ♦ Y ( J J+ 1 . k K > ♦ Y ( □ J+ 1 , K K ♦ 1 ) ) / 3 . 

flo^l 1 

lie* 


IF C J J- 1 ) 330C. 3 DOO t 3U01 

60^14 

119* 

3,00 1 

TTM2* ( TTM 1 +Y ( JJ- 1 . KK)*Y( JJ-I t‘KK*l D/3* 

0C415 

120* 


GO TO 3002 

00** I 6 

121 * 

3000 

TTM2* TTMI 

CC4 17 

1 22 * 

3002 

IF (KX-4) 30C3. 3003, 3t3o4 

00422 

123* 

3004 

TTM3* (TTM2*Y(JJ, KK* 1 ) ♦ Y ( J J ♦ 1 , KK-l))/3» 
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00*2-3 

12H* . 


. 1.£_-LJJ-*1Jl -3-00.5 *- 3.005 1 300*. 

OOH 26 

125* 

3006 

TTMM* < TTm 3*Y ( JJ-1 , KK-lD/2# 

Ot* 27 

I 26* 


GO TO 3Qin 

OCH 30 

127* 

3005 

TTMM* TTM3 

QCM31 

1 2a * 

i 

GO TO 3o 1 0 

00^32 

129* 

3003 

TTMM* TTM2 

00^33 

13C* 

. . -3 010- 

Y LJJ 1 — lUU-*-ITMM-. - . 

0C*3H 

1 3 1 * 


WRITE ( K WR , MCM9> JJ, K*. XMQi Yljj, KK> 

Q0HH2 

132* 

MQH9 

FORMAT ( 1 H • 15, 5 X , 1 5 1 2 ( 5 X 1 E12.51) 

00**3 

133* 


GO TO M 05 1 

OQHmM 

1 39* 

M0S3 

ZTEM* ABS (XNSL ( I J-ABS ( Y ( JJ| KK+nn 

00MM5 

1 35* 


If ( ZTEM-5* *XNSE<I ) ) Mq 55# M C 5 5 , MQSl 

OQMSO 

136* 

J HflSS. Z.ZJEK* lBS.l-XNS£.{.l)-ABS.CJf tJJtl . KK+IU4 

UO**S I 

1 37 • 


IF ( ZZTEM-S. *XNSE ( I ) ) M05M, M Q 5 M , M 05 1 

OCMSM 

1 38* 

M 05 1 

CONT I NUE 

aq m s a 

I 39* 

M050 

CONT I NUE 

UCHS 6 

I HO • 

C 

REWRITE CORRECTED Y ARRAY IN TEMP 

OCH40 

1 * 1 * 


IT* 1 

QGM61 

1 * 2 * 


DO- I 71 jjat. 4.4- 128 . .... . 

00^6 M 

1*3* 


Do 1 7 0 KK= 1 » 1 ?8 

0^67 

1 MM * 


TtMP ( I Y ) ■ Y ( JJ ,KK 1 

C0M70 

1*5* 


I Y * I Y + 1 

OdM7 I 

1 *6 * 

170 

CONTINUE 

OOH7 3 

1 M7* 

171 

CONT I NUE 

00*73 

I Me* 

C 

WRITE T_A.PE..RE£QBD. _ _ 

OOM75 

1 M9* 


J= 1 

OOH7 6 

1 so* 


I C NT * 1 

lid* 7 7 

IS I* 

55 

Call writer ( 9 , 2 * ier , 1 • icnt ,2om8 »temp 1 j) > 

0Q500 

152* 


IF <lCNT- 8 ) 5M, 56. 56 

aoso3 

1 5 3 * 

5 M 

J* J ♦ 20 M 8 

OdSoM 

ISM* 


ICNT*.. I_CMt . * 1- . . ._ 

00505 

1 SS* 


GO TO 55 

0Q506 

l 56 * 

56 

WRITE ( K W R » 57) IMAT 

ods 1 1 

157* 

57 

format ( 1 8 h wrote matrix no* , m 

0(S5|2 

158* 

51 

CONTINUE 

005 l H 

459* 

100 

continue 

Q05 1 6 

1 &0* 


call close 4?-*j.i... . 

2C5 l 7 

161* 


STOP 

00520 

162* 


End 


end or COMPILATION; NO DIAGNOSTICS, 

FIGURE 8 - Concluded 
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TRANSVERSE MATRIX GENERATION PROGRAM 

This is a very simple program which reads a magnetograph 
data tape, selects the U and R matrices, and calculates the transverse 
component T of the signal according to T = (U2+ R2 fz f f or each matrix 
point. The matrix of T values is then written on a new tape. The 
tape written by this program is not in the same format as the original 
tapes; it consists only of header records and T matrices. 


The data tapes for this program are the same as for the other 
programs described in this section. Two data values must be supplied 
through DATA statements; these are NMAT and DINT. The definitions 
of these variables and formats of the DATA statements are exactly 
the same as for the Noise Removal Program. Note, in particular, 
that NMAT is the total number of magnetograms on the tape, whether 
or not all those magnetograms contain transverse data matrices. On 
the output tape, however, headers and matrices are written only for 
those magnetograms which contained transverse data. 

A listing of the program is given in Figure 9. 


oomn 

I* 

00101 

2* 

00103 

3* 

00103 

H* 

001 OH 

&• 

i fin I hjl 

wu 1 Ufi 



00106 

7* 

00106 

8* 

00110 

9* 

ool ii 

10* 

00112 

1 1 • 

00 1 t~3 — — 

• 

00 I I H 

I 3* 

00 MS 

1M* 

001 16 

IS* 

00116 

16* 

00117 

17* 

— OO~1t?0 

- - 


WRITE t*pe FOR plotting TRANSVERSE field patterns 
01MENSI0N U<1638H>, RI20M8I, HT|20H8) 

DIMENSION Aj 151 I IA(RI 

NMAT* NO. OF MAGNETOGRAMS ON DATA TAPE 
DATA NM A T / 1 *4 / 

DAT-A --D J-W-T-/D.-0-A • - ■ 

D I NT *0 IF NO INTENSITY MATRIX PRESENT, .1 IF INTENSITY MATRIX 
PRESENT 
KRE« S 
KWR * 6 

CALL OPFN <8, I, M 

-C-Ai-L- REWIND ID I - 

CALL OPEN (9, 1 , b> 

CALL REWIND <91 
IMAT ■ l 

SKIP DOWN TO FIRST 0 MATRIX 

I CALL REDTPR <8,2»1ER,NW ,2 , l AM > ,7 .All I ,2*1 A<3> .iiAlfll 

• •g i ' -fO ffMA-T-f-NH •••!*•» * * 


FIGURE 9. LISTING OF TRANSVERSE MATRIX 
GENERATION PROGRAM 
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00121 19* 

00 1 28 20* 

00127 21* 

00135 22* 

00136 23* 

00M6 2*4* 

001*47 25* 

Ofrt-55 36*~ 

00156 27* 

00157 28* 

00160 29* 

00165 30* 

00166 31* 

00166 33*- 

00171 33* 

00177 3*4* 

00175 35* 

00176 36* 

00177 37* 

0000*4 38-*- 

00205 39* 

00210 <40* 

00210 *4 1* 

00210 *42* 

0021 3 *43* 

002 1 *4 4-4* 

00215 *45* 

00222 *46* 

00225 *47* 

00225 *48* 

00230 *49* 

00231 50*- 

00232 51* 

00237 52* 

002*40 53* 

002*43 5*4* 

002*46 55* 

0 03 *4 7 56*- 

00250 57* 

00255 58* 

00255 59* 

00256 60* 

00257 6 1 * 

00350 -6-3*- 

0026*4 63* 

00265 6*4* 

00265 65* 

00270 66* 

00271 67* 

00300 69* 

00301 70* 

00303 7|* 

00303 72* 

0030*4 73* 

00-305 7-4-*- 


IF ( N W - | 9 ) 1. 17, l 

\J WRITE (KWR, 21) Ntt 

26 WRITE (KWR, 167) U(l)« I A ( 2 > • All), A(2> 

167 FORMAT <1H , 2(M0, IX), 2IA12, 2* ) ) 

WRITE (KWR, 169) 1(A(I), !■ 3, 7), 1A13), 1A{4)) 

169 FORMAT (144 , 5<E12. 5,2X1, 2(110, 2X)I 
WRITE (KWR, 168) (All), 1*8, IS) 

168 FOR M AT < 1 44 i - 6-4 €l- 2- i 5-r-- ax > -) 

22 FORMAT (*4H N*«, IS, 6X,5H1CNT«, 15) 

AMG* 5,0 

WRITE (KWR, 299) 1A(2>, 1MAT, AMS 

299 FORMAT (1H , 2(|10,5X), FlO.O) 

IF ( 1 A ( 2 ) “5 ) 201, 209, 208 

X CASE 1. SK |P T O -M E 3 T H EADER 

201 CALL REDTPR < 8 , 2 , I E R , N1" , 2 , I A (1) , 7 , A ( 1 ) , 2 , 1 A ( 3 > , 8 , A ( 8 ) > 
IF ( Nim > 1 9 ) 201 , 202, 201 

202 1MAT * IMAT + 1 
AMG* 1,0 

WRITE (KWR, 299) IA(2), IMAT, AMS 
G-0- - T -0 1-7 

208 IF ( I A ( 2 ) " 6 1 ) 209, 210, 210 

209 IF (DINT) 3, 3, 203 

C CASE 2, SELECT 1ST AND 2ND MATRICES 

C CASE 3. SELECT 2N0 AND 3R0 MATRICES 

203 CALL REDTPRIB, 2, 1 ER , NW, I, ICNT, 2048, U(l>) 

AMS* 2, 0 

WRITE (KWR, 299) IA(2>, IMAT, AMS 

IF ( I C NT * 8 I 203, 3, 3 

210 IF (DINT) 203, 203, 20*4 

C CASE *4, SELECT 3R0 AND *4TH MATRICES 

204 I PE X » 0 


— A 46 * - 3 ,0 . — 

WRITE (KWR, 2991 IA(2), IMAT, AmS 

CALL 44E0TPR (8, 2, I ER , Nw , I, ICNT, 20*48 , U(l)> 

IF (JCNT-ei 204, 205, 205 

205 IF ( IOEX 1 206 , 206, 3 

206 IDEX * 1 

WRITE (KWR, 299) IA(2), IMAT, AMS 

GO TO 204 

C 3RD MATRIX on TAPE is 1 S T U MATRIX IF I Matrix IS present 


3 IU * 1 

6 CALL REoTPR (8, 2, I ER., N* , 1, ICNT, 2048. U(lU)» 


I U m IU * 2048 
IF (ICNT*8) 6, 7. 7 

HAVE STORED 1ST U MATRIX. NO* read 1ST RECORD oE R MATRIX 
CALL REDTPRl 8, 2, IER , Nw, l, ICNT, 2048, R(l)) 

WRITE (KWR, 22) N* , ICNT 


H T ( J ) * SORT ( SORT (U ( J ) *U( J ) *R (J ) *R ( J )) ) 

CONTINUE 
II, 2048 

FIRST 4 WORDS OF HT MATRIX A«E GARBAGE 
CALL WRITER (9, 2, IER, 2, 1A( | ) ,7 .Alt ) , 2 , I A ( 3 ) , 8 , A ( 8 ) ) 
XAll ILR — 2-, — |XB-t — 1-, -tCXil.,- 
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00306 75 * 

0 0 30 7 — 

00313 77 * 

00316 70* 

00317 79 * 

00320 80 * 

00322 81 * 

00^23 <^ 2 *~ — 

00320 83 * 

0032*1 84 * 

00327 85 * 

0033 ? 86 * 

00333 87 * 

0 0- 3 3 6 8 ** 

00337 89 * 

00380 90 * 

0038J 91* 

00382 92 * 


13 CALL REDTPR 18, 2, I ER , NW, 1, KNT, 2008, R<1>) 

I Sj JTE — LUJtLFL, — Z 2 i — M.ML+. . . 1. 0 1 1 

00 1 0 J* 1 , 2008 
! * J ♦ 1 1 

H T < J ) ■ SORT (SQRT(U( 1 ) *U( I ) *R( J)*R( J) ) ) 

10 CONTINUE 

II- 11 ♦ 2088 

— C AL - L *RH Tfe R — 4 T - UN T - t --^&44-,-h T < IF) 

IF CICNT-8) >3, 12, 12 

c have REACHED end of R MATRIX 

12 WRITE (K*Ri 18) IMAT 

18 format (2ih wrote ht matrix no. » i 3 1 

IF ( IMAT - NMATJ 15, 18, 18 
GO TO 1 

18 CALL CLOSE (9, 3) 

STOP 

END 


end of compilation: NO DIAGNOSTICS, 


FIGURE 9 - Concluded 
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